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ON A STOKES-TYPE SYSTEM ARISING IN FLUID VESICLE 

DYNAMICS 


DANIEL LENGELER 


Abstract. This article is the first in a series of papers on the analysis of a 
basic model for fluid vesicle dynamics. There are two variants of this model, 
a parabolic one decribing purely relaxational dynamics and a non-parabolic 
one containing the full dynamics. At the heart of both variants lies a linear 
elliptic system of Stokes-type. Understanding the mapping properties of this 
Stokes-type system is crucial for all further analysis. In this article we give 
a basic exposition of the dynamical model and a thorough L 2 -analysis of the 
Stokes-type system that takes into account geometric variations of the fluid 
vesicle. 


1. Introduction 


The basic constituent of most biological membranes is a two-layered sheet of 
phospholipid molecules, a lipid bilayer. Phospholipid molecules possess hydrophilic 
heads and hydrophobic tails. When exposed to water they arrange themselves into 
a two-layered sheet with tails pointing inward. Due to the hydrophobic effect, 
biological membranes tend to avoid open edges and form closed configurations, 
the vesicles. While a membrane is only a few nanometers thick, the diameter 
of the vesicles formed can be up to 10^ times larger. This separation of length 
scales suggests to describe vesicles as two-dimensional surfaces embedded in three- 
dimensional space. Since the solubility of the phospholipids is very low, the ex¬ 
change of material between the membrane and the ambient solution is negligible. 
Thus, vesicle configurations are not determined by a surface tension but rather by a 
bending elasticity leading to an amazing variety of equilibrium shapes; in contrast 
to the typically spherical equilibrium shapes of fluid interfaces. A basic model for 
the elastic energy, the Canham-H elf rich energy 
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was proposed independently in three seminal papers [smiiiii]. Here, F is the 
two-dimensional surface representing the membrane, H and K denote twice its 
mean curvature and its Gauss curvature, respectively, and k, kq are the bending 
rigidities with the physical dimension of an energy. Furthermore, the constant Cq 
is the spontaneous curvature which is supposed to reflect a chemical asymmetry 
of the membrane or its environment. The Canham-Helfrich energy (with Cq = 0) 
can be seen as the lowest order term in an expansion in curvature. By the Gauss- 
Bonnet theorem, the Gaussian part of the energy is a topological constant and 
hence irrelevant for vesicles of fixed topology. Due to the above mentioned low 
solubility of the phospholipids and a low compressibility of the membrane its area 
is practically fixed. Furthermore, despite the fact that the membrane is permeable 
to water, the fluid volume enclosed by the vesicle is also fixed, due to strong osmotic 
effects. 

Thus, in the simplest case, equilibrium vesicle shapes supposedly minimize the 
Ganham-Helfrich energy under the constraints of fixed area and enclosed volume. 
The known solutions to this minimization problem can indeed be observed exper¬ 
imentally. However, other (non-axisymmetric) shapes found in the laboratory like 
the starfish vesicles and shape transformations like the budding transition seem 
not to be contained in this model. The reason for this deficiency is that the bilayer 
architecture of the membrane is not correctly incorporated. At a characteristic 
temperature all lipids transition from a gel to a liquid phase. The lipid bilayers we 
are dealing with are assumed to be in the liquid phase, allowing the monolayers 
to freely flow laterally and to slip over each other while the membrane retains its 
transverse structure. This leads to a lateral adjustment of the monolayer densities 
to curvature. Incorporating the bilayer architecture in the form of inhomogeneous 
monolayer densities yields, to lowest order, the additional elastic energy 

(2) FsL = ^J ((P+ - dH)^ + (p- + dHf) dA, 

r 

introduced in [39) . Here, are the reduced density deviations of the monolayers 
and d is the membrane thickness. By scale considerations one can show that in 
equilibrium the mean density {p^ + p~)/2 is a fixed constant leading again to 
a fixed membrane area while the density difference {p'^ — p~)/2 gives a nonlocal 
energy involving the total mean curvature. Together with the volume constraint this 
yields the area-difference elasticity model introduced independently in gniiiiiE]. 
This model seems to be quantitatively consistent with all experiments performed so 
far. An intermediate model, the bilayer couple model, introduced in [njiiiiis], is 
based on the Canham-Helfrich energy, but takes into account three hard constraints, 
namely area, enclosed volume, and total mean curvature. 

In most experimental situations the dynamics of fluid vesicles are of interest, 
too. On the one hand, thermal excitation always leads to fluctuations around stable 
equilibrium shapes. On the other hand, real biological membranes are constantly 
put into non-equilibrium states by biophysical processes, like molecule insertions 
and extractions, or by external hydrodynamic flows. When tweaking a parameter 
(like the enclosed volume, e.g. by changing the osmotic conditions) discontinuous 
shape transitions may occur, that is, a formerly stable shape becomes unstable and 
decays towards a new minimum. In any case, it is crucial to take into account 
the interaction of the moving membrane with the viscous bulk fluid. Apart from 
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the bulk viscosity, two other relevant dissipation mechanisms have been identified, 
namely friction between the two monolayers and shear viscosity within each layer; 
see for instance mM- In most situations of interest both the bulk and the surface 
Reynolds numbers are very small, typically of the order In such situ¬ 

ations it makes sense to neglect inertial effects in the bulk and on the membrane, 
leading to purely relaxational dynamics. The significance of the bulk fluid then lies 
in the induced nonlocal self-interaction of the membrane and in the viscous energy 
dissipation. 

Several models for the dynamics of fluid vesicles have been proposed; see for 
instance [34l[39l[T5l[28l[24l[2l[32]. Some of these are restricted to the minimal energy 
o while others take into account the bilayer architecture using the additional 
energy (ED. While |39] deals only with an almost planar membrane, [MIIIIIEHI 
124] neglect in-plane shear viscosity. Including the latter in an arbitrary geometry 
leads to the well-known Boussinesq-Scriven surface fluid [ZIISS] as basic continuum 
mechanical modeling shows; cf. [2132]. In this paper, we take into account the 
in-plane shear viscosity while neglecting the bilayer architecture of the membrane. 
More precisley, we study a single homogeneous Newtonian surface fluid subject 
to additional stresses induced by the Canham-Helfrich energy, interacting with a 
homogeneous Newtonian bulk fluid. The full model (and intermediate models) will 
be studied in future publications. For a comprehensive presentation of the physics 
of fluid vesicles we refer the reader to m and the references therein. 

Phase diagrams of the different curvature models assign the shape of lowest 
energy to each choice of external parameters (like the enclosed volume). Most 
of the insight into these phase diagrams and into the bifurcation and stability 
of stationary shapes has been generated by numerical computations and formal 
(linear) analysis; again, see |37j and the references therein. From the point of view 
of a rigorous mathematical analysis, not much is known concerning these questions. 
Critical points of the energy CD subject to the volume and area constraint solve 
the Helfrich equation 

(3) grad^2 F’ -I- Ai -|- X 2 H — 0 

with Lagrange parameters Ai,A 2 . While the round sphere is a solution of this 
equation for any choice of Cq, it is the only known analytical solution of spherical 
topology for Co = 0; see [23] for further analytical solutions of spherical topology 
in the case of non-vanishing spontaneous curvature. In [29] the authors prove the 
existence and analyze the stability properties of a one-parameter family of critical 
points bifurcating from the round sphere. In |35] it is shown that there exists a 
global minimizer of the Canham-Helfrich energy with vanishing spontaneous curva¬ 
ture in the class of smoothly embedded topological spheres with fixed isoperimetric 
ratio a € (0, i]Q Furthermore, it is shown that the minimal energy is a continuous 
and strictly decreasing function of cr, approaching Stt from above as cr \ 0 and 
that, roughly speaking, the minimizers converge to a double sphere in a measure- 
theoretic sense in this limit. The existence result is generalized in m to higher 
genus surfaces. Experiments indicate that all minimizers are axi-symmetric (even 
for non-vanishing spontaneous curvature); however, proving this mathematically 
seems to be a hard task. Concerning the more realistic models, nothing is known 


^Since the Canham-Helfrich energy is scaling invariant, instead of fixing both the enclosed 
volume and the area, it is sufficient to fix only one parameter, for instance the isoperimetric ratio. 
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from the point of view of a rigorous mathematical analysis. Insight into the dy¬ 
namics of fluid vesicles based on numerical computations and formal linearized 
analysis can be found, for instance, in [531 IMl [HI UHl I3H1 IMl 1211311321131131 IS] ; h 

seems that [4] is the first article presenting fully three-dimensional numerical com¬ 
putations for the model we are studying in the present article. However, again, 
not much rigorous analysis has been done on this topic. Concerning the Canham- 
Helfrich flow, that is, the L 2 -gradient flow of the Canham-Helfrich energy with 
prescribed enclosed volume and area, a partial local well-posedness result has been 
shown in |5^. There exist further results [UlllllTl concerning a Helfrich-type 
flow where the Lagrange parameters instead of volume and area are prescribed and 
which consequently should not be related directly to fluid vesicles. In m local-in- 
time existence and uniqueness for a homogeneous Newtonian surface fluid subject 
to Canham-Helfrich stresses is shown. While the bulk fluid is neglected the au¬ 
thors keep the inertial term in the equations for the surface fluid, yielding a kind 
of dissipative fourth order wave-type equation. In m local-in-time existence and 
uniqueness of a homogeneous Newtonian bulk fluid with inertial term interacting 
with a compressible, inviscid surface fluid without inertial term is shown in the 
L 2 -scale, the membrane model being rather non-standard. 

The present article is the first in a series of papers on the relaxational and 
oscillatory (that is, involving both inertial terms) dynamics of a homogeneous sur¬ 
face fluid subject to Canham-Helfrich stresses and interacting with a homogeneous 
bulk fluid. It contains a basic exposition of the dynamical models and a thorough 
L 2 -analysis of a linear Stokes-type system that plays a fundamental role in these 
models. Understanding the mapping properties of this system is an essential first 
step for all further work on topics like local well-posedness, stability, and long-time 
behavior. We will see that our Stokes-type system induces a Riemannian metric 
on the manifold of surfaces (of fixed area and enclosed volume) that is part of a 
gradient flow structure for the purely relaxational dynamics; this observation will 
prove very useful for stability analysis. For most further investigations, in fact, it is 
important to establish an L 2 -theory that holds uniformly with respect to suitably 
controlled membrane deformations. This is particularly important and delicate 
for the analysis of the formation of singularities. It is a well-known and obvious 
fact that the Canham-Helfrich energy is invariant under isotropic scaling which is 
a continuous, non-compact symmetry. Although this symmetry is broken by the 
constraints of fixed area and enclosed volume, one can expect that a blow-up of 
curvature can only be ruled out by proving or assuming non-concentration of en¬ 
ergy; see for instance [44] for the general heuristics on energy-critical systems]^ A 
similar result has been proven for the Willmore flow in [21] . The proof of such a 
fact relies on higher order energy estimates; hence the need for an L 2 -theory of the 
underlying Stokes-type system that involves a precise control of the dependence on 
membrane deformations. 

The present article is organized as follows. In Section[2]we shall derive the math¬ 
ematical models and perform some basic computations; in particular, we present the 
gradient flow structure underlying the purely relaxational dynamics. In Section |3| 

^Here, additionally, we need to assume or prove (which would be an interesting task) that no 
self-contact of the membrane occurs. Furthermore, the Canham-Helfrich energy, involving only 
the mean curvature, lacks coercivity. Thus, in order to rule out blow-up of curvature, we will have 
to control the integral of the second fundamental form squared instead of the mean curvature 
squared. 
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we perform the L 2 -analysis of the Stokes-type system for smooth, fixed membranes 
and then for variable membranes. In Appendix |A] we recall some basic facts con¬ 
cerning the kinematics and dynamics of fluid interfaces, and we derive the stresses 
resulting from the Canham-Helfrich energy. In Appendix [B] we present some useful 
formulae concerning covariant differentian and curvature, while in Appendix [Cl we 
prove some facts on the function spaces we shall use. 

Before we proceed, let us fix some notation. Throughout the article (apart from 
the appendix), we denote by e the Euclidean Riemannian metric in K^, and for 
surfaces E C we denote by g the metric on T induced by e. We also use 
the notation u ■ v instead of {u,v)e for u,v G We denote by Pr the held of 
orthogonal projections onto the tangent spaces of T and by TT the tangent bundle of 
r. Furthermore, [ujp denotes the trace of the bulk held it on T; however, when there 
is no danger of confusion we will sometimes omit the brackets. Moreover, we write 
k, H, and K for the second fundamental form, twice the mean curvature, and the 
Gauss curvature of T with respect to e, respectively. With a slight abuse of notation 
we use same symbol k also to denote the Weingarten map, that is, in coordinates 
we write kap and fcf. Furthermore, for any metric e on an arbitrary manifold, we 
write V®, gradg, fcg, etc. for the associated Christoffel symbols, differential 
operators, and curvature terms, and we use the abbreviations F^^ := V := V®, 
grad := gradg, etc. for the corresponding Euclidean objects. When working in 
coordinates and confusion about the underlying metric can be ruled out, we use 
the semicolon to separate indices coming from covariant differentiation from the 
original indices; for instance, for a covector held w we write = uJi-j. Again 

with a slight abuse of notation, we write (Ti,T 2 )g for the contraction of tensor 
helds Ti of the same rank using the metric e and |T|e for the corresponding norm. 
For manifolds N,M, a. diffeomorphism (p : N ^ M, and tensor helds T on N and 
S' on M we denote by the pushforward of T and by (p*S the pullback of S 
under p. We denote by r(a) generic tensor helds that are polynomial or analytic 
functions of their argument a such that r(0) = 0. Furthermore, for tensor helds 
ri and r 2 we write ri * r 2 for any tensor held that depends in a bilinear way on 
ri and r 2 , and we use the abbreviations r * {ri,... ,rk) = r*ri + ... + r*rk and 
= r * ... * r (with k factors on the right hand side). For s > 0, p,q G [1, oo], and 
tensor helds T we denote by ||T||p, \T\s^p, and \T\s^p^q the corresponding Lebesgue, 
Sobolev-Slobodetskij, and Besov (semi-)norms, see Appendix ICl we omit the spatial 
domain in the notation since it will always be identical to the domain of dehnition 
of the argument T. We denote by and the L 2 -scale of Sobolev-Slobodetskij 
spaces, where the lower index 0 refers to vanishing traces (for s > 1/2). 


2. Modeling and basic observations 


2.1. The mathematical models. Let Id be a bounded domain in containing 
a homogeneous Newtonian huid and a closed lipid bilayer Ft depending on time t. 
We have the usual Navier-Stokes system in the bulk il\Tt'. 


Du 

rt.^ = d„S, 

div u = 0. 


Here, u is the huid velocity, = dt ■ -l-(u • V) • is the material derivative, in 
particular, ^ is the huid particle acceleration, S = 2g,i,Du — T:I is the (bulk) stress 
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tensor, Du is the symmetric part of the gradient of u, tt is the pressure, pb is the 
constant bulk fluid density, and pt is the constant dynamic viscosity of the bulk 
fluid. Obviously, we have div^" = pb^u — gradvr. We assume that u vanishes on 
dfl. As is done in all dynamical models mentioned above, we prescribe the no-slip 
condition at the membrane, that is, the velocity of the bulk fluid at the membrane 
equals the membrane velocity. For the normal part of the velocity, this can be 
justified by assuming an instantaneous osmotic equilibrium of the bulk fluids on 
either side of the membrane. Mathematically, the no-slip condition implies that 
the fluid velocity is continuous across the membrane. On Ft we decompose the 
membrane velocity u into its tangential and normal part, u = v + w v, where v is 
the outer unit normal. In the following, we have to consider hybrid tensor fields 
on Ft, that is, tensor fields that involve both tangential and non-tangential vectors 
and covectors; see for example [Mill]. An instance of such a hybrid tensor field is 
the velocity field u on Ft that is non-tangential in general. Another instance is the 
surface stress tensor T which is a hybrid (1, l)-tensor, taking a tangential direction 
and returning a force density that is, in general, not a tangent vector. We adopt 
the standard coordinate notation for such tensors and write T^, the greek index 
taking the values 1,2 and the latin index taking the values 1,2,3. In this paper, 
the latin index will always refer to Cartesian coordinates in while the greek 
index refers to arbitrary coordinates on F^. Likewise, we write u® for the Cartesian 
components of (not neccessarily tangential) vector fields u on Ft. There exists a 
canonical covariant differential calculus for hybrid tensors; see [1]. Here, however, 
we only need the divergence operators, which take the form 


Divu = g°‘^{daU,di3)e, 

{mvTy = g‘^^Tf,.^p, 


where the semicolon in the second line denotes the usual covariant differentiation 
of the covectors (T^)a=i ,2 (for fixed i). Basic continuum mechanical considerations 
for a general interfacial fluid, see Appendix or [36l [T], lead to the following 
conservation laws for linear momentum and mass on Fj: 



( 4 ) 


Here, p is the surface fluid density and [S'] is the jump of the bulk stress tensor 
across the membrane (subtracting the outer limit from the inner limit). We assume 
p to be constant, reducing (l4])2 to Divit = 0. Furthermore, assuming the interfacial 
fluid to be isotropic and Newtonian, the (strictly tangential) fluid part 
of the surface stress tensor takes the form 



see m^- Here, q is the surface pressure, acting as a Lagrange multiplier with 
respect to the constraint Div u = 0, p are the surface dilatational and shear 
viscosities, respectively, and Du is the surface rate-of-strain tensor. Analogously to 
the bulk case, Du is the Lie derivative of g of with respect to u, {Du)ap = ^{va-,/3 + 
yp;a) — u) kap; SCO Appendix lAl or [Mill]. We have g'^^{Du)-yb = divg u — wH = 
Divu = 0; for the second identity see below. We conclude that the fluid part of the 
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surface stress tensor takes the form 

= -qS^+ 2fl [VuYa = {v a-^ + v^. a “ 2wka^). 

The Canham-Helfrich energy m induces stresses in the membrane which are given 

by 

"“T; = k{{H- C^f/2 d\-(H- Co) kid} -{H- Co),aC), 

where the comma denotes the usual partial differentiation; see Appendix |A] or [lOj . 
We have 

(5) = -K{AgH + H{H^/2-2K) + Co(2K - HCo/2))iy = -gmdL^Fiy. 

Finally, our constitutive assumption concerning the interfacial fluid is 


The full system of equations now reads 


w|^ = divS 

in D \ Ft 

div u = 0 

in D \ Ft 

p^ = Divr + [5iu 

on Ft, 

Div u = 0 

on Ft, 

u = 0 

on dfl. 


This is the non-parabolic model we shall study. It can be interpreted as a vector 
field on the nonlinear manifold consisting of embedded surfaces T C of hxed 
area and enclosed volume and of velocities u : —>■ in the sense that (r,u) 

is mapped to ([ujr • v,dtu)^ where dtu is computed from (ISl)i, 3 ; analogously to the 
classical Navier-Stokes system one can show that the pressure functions tt and q at 
some point in time can be computed from the knowledge of T and u at that same 
instant alone. 

Let us put the terms involving the surface divergences in a more explicit form. 
From the decomposition u = v + w v we infer 

= d^v‘^ +v^ST'^l^-wki = dWgV-wH. 

Furthermore, decomposing the surface stress tensor into a tangential and a normal 
part, = ^T^d} + "^TaC, we obtain 

(DivT)* = 

Note that d} is the a-th component of the covector X H> {X, e*)e, where e* is the i- 
th standard basis vector in Thus, d}.^ = d} „ = kpaC. Furthermore, 

by definition of the Weingarten map, we have C= —kid}. Hence, we obtain the 
general formula 

(7) (Div ry = (divg ‘T)“5; + *T"^ko,pC + (divg "T) C - ^T‘^kid}. 

From this, let us compute Div ^T. We have 

(divg^T)^ = -g°'^q,a + g.g^^g^'*{va-js + v^.^aS - 2{wka^).^s)- 
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By definition of the Riemannian curvature tensor R, we have Va-,-yS = Va^Sj + 
Recalling that for a surface Rg^^^ = K {S^g^s — S^ga-y) we obtain 
g°‘^Va-,~iS = (divgu),^ + Kvry. Smco div^u = w H, this gives 

divg = — gradg q + g (AgU + gradg(?« H) + Kv — 2divg(r(; A:)). 

Furthermore, we have 

= -qH + 2gi {Vu)^^k^p 
= —q H + 2fj, {ya-pk°'^ — w {H^ — 2K)). 

In view of © and (O we conclude that 

Div T = — gradg q — q Hv + g (AgU + gradg {w H) + Kv — 2 divg(?ii k)) 

(8) +2g{{V<^v,k)g-w{H^ -2K))iy 

- K{AgH + H{H^/2 - 2K) + Co{2K - HCqI2))v. 

Note that the highest order term in the normal part of Div T is a quasilinear fourth 
order differential operator acting on the membrane configuration while the highest 
order term in the tangential part is a second order operator acting linearly on the 
tangential membrane velocity. Recalling that Du/Dt is the fluid particle acceler¬ 
ation, we see that ®3 is a fourth order wave-type equation governing membrane 
deformation which is strongly coupled to a two-dimensional Navier-Stokes type 
system governing the lateral flow on the membraneH 

Let us have a look at the order of magnitude of the bulk Reynolds number TZi, and 
the surface Reynolds number TZ of the system ® in typical experiments; cf. [37) . 
Typical relevant length scales are of the order Ltyp = Igm, while typical time scales 
(accessible to video-microscopy) are of the order Ttyp = 10“^s. Furthermore, typical 
values for the bulk density and viscosity are pt = lO^kg/m^ and pb = 10~^kg/ms, 
respectively. This gives 

Ub = = 10-3. 

P-bTtyp 

Assuming that the membrane volume density is comparable to that of water, for 
its surface density we obtain p = lQ~^kg/m?. Typically, the surface viscosity is of 
the order g = \Q~^kg/s. Thus 


nT'^ 

= aAAp = 10-5. 

g Ttyp 

Hence, under such circumstances we can expect that the viscous forces strongly 
dominate the inertial forces. Neglecting the inertial terms in ®, we obtain a set 


^More precisely, the equation governing membrane deformation is not hyperbolic, but it is 
dispersive. The principal part of its linearization has the form p9j-|- kA^. This operator essentially 
describes the linear evolution of a thin elastic plate. Note that it factorizes into two Schrodinger 
operators; in particular we can expect an infinite speed of propagation of disturbances. 





ON A STOKES-TYPE SYSTEM ARISING IN FLUID VESICLE DYNAMICS 


9 


of equations describing purely relaxational fluid vesicle dynamics: 

div S = 0 in fl \ r t , 

div u = 0 in fl \ r t, 

(9) Divr = -[5'|u onTt, 

Div u = 0 on Ft, 

u = 0 on (9D. 

Here, phase space consists of the embedded surfaces F C D of fixed area and 
enclosed volume alone. Note that, by ([5]), equation @3 can be written in the form 

(10) Div = grad^^ F v. 

Thus, at a fixed instant in time we have to solve for the fluid velocity u a linear 
Stokes-type system with prescribed Neumann-type boundary values given by the 
L 2 -gradient of the Canham-Helfrich energy. Mapping F to \v\r ■ v then defines 
the dynamics of our system. Compared to the classical Canham-Helfrich flow (or 
Willmore flow), there is an additional Neumann-to-Dirichlet operator involved here. 
Since the L 2 -gradient of F" is a fourth order operator, the mapping F 1 —[it]r • v can 
be considered as a nonlinear, nonlocal pseudo-differential operator of third orderQ 
The Stokes-type operator defined by the left hand sides of (l 9 ])i, 2 , 4,5 and ITOl) will 
also play a crucial role in the analysis of ([ 6 ]). 


2.2. Some basic observations. Let F), z = 1,..., fc, denote the connected com¬ 
ponents of Ft, and denote by D) the part of D enclosed by FJ. Reynolds’ transport 
theorem then gives 




= J div udx = 0. 
oj 


In particular, the volume enclosed by each connected component of the membrane 
is conserved. Similarly, we can show that the area of each connected component is 
conserved: 





vdA = 0. 


Pt 


Now, let us present the gradient flow structure underlying the parabolic flow ®. 
Consider the manifold N of embedded surfaces of fixed area and enclosed volume. 
For T G N, the tangent space TpN can be identified with the space of scalar fields 
w on T such that the linearized constraints 


( 11 ) 


? dA = 0 and 


w H dA = 0 


pi 


Y'i 


“^The situation is analogous to the relation between volume-conserving mean curvature flow and 
the Mullins-Sekerka system. There, however, we have a Dirichlet-to-Neumann operator involved. 
The mean curvature operator being of second order, the Mullins-Sekerka gradient is a third order 
operator, as in the present case. 
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hold for all connected components F® of F. For w S TpiV, consider the system 

div S = 0 in n \ F, 

div u = 0 in n \ F, 

Pr(Div/T+[51i.) =0 on F, 

Div u = 0 on F, 

u - V = w on F, 

u = 0 on do,. 


( 12 ) 


Note that the conditions dm are necessary for the solvability of these equations, 
due to the incompressibility constraints. For wi,W 2 G TrN, define the Riemannian 
metric on N associated with fluid vesicle dynamics by 

(13) (ici, 102 ) E := 2/16 J {Dui,Du2)edx + 2fj, J {'Dui,'Du2)g dA, 

n\r r 

where Ui , U 2 solve the system m with data Wi,W 2 - Note that the length of a 
curve in N endowed with this metric is given by the energy dissipated during the 
corresponding forced deformation of the membrane. The representation of —dF 
with respect to the metric (IT^ is given by [u]r • i', where u solves dH). Indeed, for 
all w GTrN and corresponding solutions u of dD we have 


{u ■ i',w)v = J {Du, Du)f. dx + 2fi J ('Du,'Du)gdA 
o\r r 

( 14 ) = - J {lSjv,u)edA - J {Div fT,u)edA = J {Div^T,u)edA 

r r r 

= —J F w dA =—dF{w). 

r 

Here, S and denote the stress tensors with respect to u, and we used integration 
by parts for the second identity (see below), dUa for the third identity, and di]) for 
fourth identity. We conclude that, indeed, m is the gradient flow of the Canham- 
Helfrich energy on N endowed with the Riemannian metric m- In particular, the 
energy F is a strict Lyapunov functional, and, along the flow, 

^F = dF{u ■ v) = —{u ■ v,u- v)v = — 2/16 J \Du\l dx — 2^ J \Vu\^ dA. 

n\rt Et 

It remains to prove the second identity in (ITT)) . In view of dHl)i ,2 we have 

0= f {div S,u)edx = — j {\S\i', u) e dA — f {S,\7u)edx 
n\r r n\r 

= — J{lSli',u)edA —2fj,b J {Du,Du)edx 
r n\r 

Furthermore, we compute 

= -qg^P{d^,dpu), + 2gg<^f^{Vu)l {d^,dpu)e. 
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Note that g°'^{da,di3u)e = Divu = 0 and {d-y,dpu)e = vp-a — wkap- Using inte¬ 
gration by parts and the fact that "Du is symmetric, we obtain 

J (Div^T, u)edA = —2g, j{'Du,'Du)g dA. 
r r 

Furthermore, from the computations above and Reynolds’ transport theorem we 
see that the energy identity for ([6]) reads 

+f y + j \Du\ldx-2g, j\V u\ldA. 

n\rt Tt n\rt Ut 

Finally, let us analyze the equilibria of ([5]) and (jH]) a little more closely. From 
the energy identities we see that Du = 0 in U \ F; taking into account the no-slip 
conditions and Korn’s equality, this gives u = 0 in U. Now, m 1 and m 1 reduce 
to gradTT = 0 in fl \ F, that is, tt is constant in each connected bulk component. 
Furthermore, (HJa and ©a reduce to 

gradg q + qHv -b v = -|7r]i/. 

Since the surface pressure gradient is the only tangential component in this equa¬ 
tion, q is constant on each connected component of F. Thus, we have 

gradi^ F-b [tt]- b g i/= 0. 

This is the Helfrich equation ([3]) with the pressure jump and the surface pressure 
acting as Lagrange multipliers with respect to the volume and area constraints. 

3. L 2 -THEORY FOR Stokes-type systems 

Let U C be a smooth, bounded domain and F C U a smoothly embedded 
closed surface. We write F*, i = 1,..., fc, for the connected components of F, 
for the open set enclosed by F*, and we let 

k 

■.= n \ (|Jf*uu*). 

i=l 

Let us assume that /ib,/r > 0. We shall analyze two slightly different Stokes-type 
systems. The first one is related to the Neumann-to-Dirichlet operator discussed 
above and reads 



div S = fi 

in D \ F, 


div u = f 2 

in U \ F, 

(15) 

Div/r + |,s]p = /3 

on F, 


Div u = fi 

on F, 


u = 0 

on dQ 

for appropriate data (/i. 

..., fi). From (fT^9 

we obtain 
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Thus, recalling Div u = divg v — w H, by (llhll /i we have 

(16) J fi/H dA = — J f 2 dx for each T* that is a CMC surface, 
r* n* 

Recall that the only closed CMC (= constant mean curvature) surfaces embedded 
in are (collections of) round spheres. Furthermore, of course, we have 


(17) 


/2 dx = 0. 


Hence, for s > 1, we define the space of data 

F«(r) := {(/i,..., h) I h e \ r;K"), /2 e \ r), 

Pvh e i7'*-'(r;Tr), u e 

such that (fTBll and (fT7)l hold} 

endowed with the canonical norm. Furthermore, for s > 0, we introduce the space 
of solutions 

]E"(F) := {(M,7r,q)|uG £C+i(H\r;R3)n77i(fl;R3)^ p^[u]r G FC+i(F;TF), 

71 G H^{Q \ F), g G i7®(F) such that (fTBl) . (flTl) hold with 

/2 = TT and /4 = q} 

endowed with the canonical norm. The conditions on tt and q provide a gauge 
fixing; as is typical for Stokes-type equations the pressure functions in (1151) are not 
uniquely determined. The second Stokes-type system we shall analyze is related to 
the Riemannian metric discussed above and reads 


(18) 


div S = fi 

in 0 \ F, 

div u = f 2 

in 0 \ F, 

PriDivfT+lSP) = f3 

on F, 

Div u = f 4 

on F, 

u-v = h 

on F, 

u = 0 

on 90 

Of course, again, (fTH)) and (fT7|) must hold. 


(19) J f 2 dx = J f^dA and J /4 -I-/s dH = 0 for alH = 1,..., fc. 

Qi pi p" 

Note that m is contained in (HU). Hence, for s > 1, we define the space of data 
:={(/i,...,/5)|/ieH«-i(fl\F;K3), /2 G \ F), /a G i7«-i(F; TF), 

/4 G i7"(F), /s G 

endowed with the canonical norm. Furthermore, for s > 0, we introduce the space 
of solutions 

E: (F) := { {u, TT,q)\ue \ F; K3) n (fl; R^), Pr Mp G (F; TF), 

TT G i7®(H \ F), g G i7®(F) such that (fT91) holds with 

/2 = tt, /4 = g, and /s = 0} 

endowed with the canonical norm. 
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For the applications we have in mind it is important to include the case of s 
being half-integer valued. For instance, both the local well-posedness of ([6|) and 
the analysis of the formation of singularities in © rely on energy estimates on the 
surface level, that is, yielding an integer class regularity for the normal velocity 
[it]r ■ i'- But then we must apply the theory developed in the present article with 
s = n + 1/2 for some n G N. Since it doesn’t cause any additional effort we prove 
our results for arbitrary real s > 1. 

Throughout the article (except for a few formulae in the proof of Theorem 13.131 
and in Appendix [C|) , we keep track of the physical dimensions in all mathematical 
expressions. This mainly serves as a rough test for the validity of the estimates we 
prove. We write ci, I G h, for generic positive constants with the physical dimension 
of length* depending only on 17, F, the Saffman-Delbriick length fi/fib, and s; let 
c := Co. 

3.1. Analysis for a fixed membrane. Let us start with the analysis of (fT51) . To 
begin with, let us assume that /2 and vanish, and define 

X := {uG I divu = 0 in 17 \ F, Divu = 0 on F, PrMr £ H^T-TT)}. 

Multiplying (TTSt i by yi S A, integrating over 17 \ F, and integrating by parts we 
obtain 


( 20 ) 


0 = J (div S, (p)e dx - y (/i, V5)e dx 
Q\r n\r 

= - f {S,X(p)edx - f {lSliy,ip)edA- f {fi,ip)edx 
n\r r n\r 

= - j {S,Dip)edx + J(DiyfT,ip)edA-J{f3,ip)edA-J {fi,ip)f,dx 


n\r 


n\r 


= —2fj,b J {Du, Dip)e dx — 2fi J {'Du,'Dip)g dA 
n\r r 

- j{f3,^)edA- J {fi,ip)edx 

r n\r 

=: -B{u,ip) + F{ip). 

Let us have a closer look at the bilinear form B. 

Lemma 3.1. B is a symmetric, continuous, and coercive bilinear form on X. For 
all u G X, we have 


( 21 ) 


A||u||^ < l|Vu ||2 < —B{u,u), 
||V®u||i< -P(u,u), 


where A is the smallest eigenvalue of the Dirichlet-Laplacian in 17. 
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Proof. Symmetry and continuity are obvious. Let us prove coercivity. Recall that 
2\Du\l = |Vu|g + div((u • V)u). Thus, we have 


( 22 ) 


J \Du\1dx = 2 j \Du\^^dx = J \yu\1dx. 


Q\r 


By Poincare’s inequality, m 1 follows. Moreover, we have 

\Vu\l = I Def(u)|g + w'^\k\l - 2(Def(i;), k)g w, 

where T)ei{v)ap = \{va-,p + vp-a) is the deformation tensor of v, and, similarly to 
the bulk case, 2| Def(u)|g = |Vu|g + divg(V®u) — However, by definition 

of (intrinsic) curvature, we have Vp.^v°‘ = (divi;),o,r’“ + Thus, taking into 

account div^ v = w H, we obtain 


2 J \Vu\ldA = J {\Vv\l+w^(H^ + 2\k\l)-K\v\l-A{Dei{v),k)gw) dA. 
r r 


This gives 

J Ni’\ldA<-B{u,u)+ \\k\\l^ J \v\ldA + 4:\\k\\^ J \Vv\g\w\dA. 
r r r 


Applying Holder’s and Young’s inequality, the continuity of the trace operator, and 
Poincare’s inequality, we obtain 


(23) 


J iVul^dA < ^H(u,u) + 16||fc|l^ J \v\l+w^ dA 

r r 


< —B{u, u) + c_i [ 

M J 


|Vu|g dx. 


Combining (l^ i and (1^ we obtain (1^ 9. 


□ 


By Riesz’ representation theorem in Hilbert space, for each F G X' there exists 
a unique u G X such that B{u, (p) = F{ip) for all (p G X. Now, we reconstruct the 
pressure functions. Consider the space 

Y-.= {uG Hl{n- R3) I p^[u]r G H\T- TP)} 
endowed with the norm 

||u||^:= ||Vu||2 + |H|^||V®PrMr||^ 


and the space 

Z := {(/2,/4) G L 2 {^) X L 2 (r) such that IH]), ([T7]) hold} 
endowed with the norm 

ll(/2,/4)|||:= ||/2||^ + |f^|^||/4|ir 

Note that X is a closed subspace of Y. Each (tt, q) G Z defines an element of Z' 
via 

(/ 2 , J f 2 dx + J q U dA. 

n\r r 
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In fact, by Riesz’ theorem, this defines an isometry 

(24) 

where 

lll(^,g)lllz:= + 

defines an equivalent norm on Z. 


Lemma 3.2. The oyerator div : Y/X —>■ Z, defined by 

ip !-)■ (div(/j, Div[(/?]r), 

is an isomorphism. In particular, we have 

(25) \\p\\y/x <c\\T)\vlp\\z. 


Proof. Obviously, div is well-defined and injective, so let us prove surjectivity. To 
this end, let (/ 2 ,/ 4 ) S Z. Choose a function w on T in the following way: On T® 
choose w G span{l,iJ} such that 


(26) 


J wH dA = — J 

pi pi 


/4 dA and w dA = 1/2 dA, 


which is possible by definition of Z. Now, let v G H^{T; TT) solve 


divg V = f 4 + wH 

on r®; note that the right hand side is mean value free and employ, for instance, 
L 2 -theory for the Laplacian on T®. Let ip := v + w G and extend this 

function to the bulk, giving po G iLo(0;K^). Furthermore, by L 2 -theory for the 
divergence operator in O \ F there exists a. pi G FIq (O \ F) such that 

(27) div Pi = f 2 - div po 

in fl \ F; note that the right hand side is mean value free in O® for alH = 0,..., fc. 
For i 0 this is obvious from (1^511 while for i = 0 we additionally take into account 
(HZl) and compute 


J div Po dx = — J w dA = — J /2 dA = J f 2 dA. 
no r UiD* no 

Finally, for p := po + pi G Y we have div ly = (/ 2 , fi). 


□ 


Corollary 3.3. The operator V : Z — C Y', defined by 

V (tt, q) (p) = — J TT div pdx — J q Div p dA 
n\r r 

for p gY , is an isomorphism. We have 
(28) lll(7’',g)|||2 < c||V(7r,g)||r. 

with the same dimensionless constant c as in 


Proof. This follow from Lemma 13.21 the fact that 

V = -div' : Z' {Y/X)' 

is an isomorphism, from the isometry (Y/X)' ~ X-^, and by Riesz’ isometry 
Recall also that adjoining preserves the operator norm. □ 
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Note that Corollary 13.31 and Lemma 13^ are equivalent to the inequality 
inf sup --^ 


IE 


^ y TT div + J q Div ip dA^ > 0, 


n\r 


which is a Ladyzenskaja-Babuska-Brezzi-type condition', cf. [1]. 

Consider the subspace U C L 2 {k}) x L 2 {T) defined as follows: {n,q) S 17 iff for 
all i = 1,..., fc we have 

(i) TT = Ni in n®, TT = kq in 9 = N® on L® with Ki, kq, S K. 

(ii) If r® is a round sphere with H denoting twice the mean curvature, then 
Ki — kq = hdH. 

(iii) If r® is not a round sphere, then k® = 0 and Hi = /Co¬ 
lt is not hard to prove that {L 2 {k}) x L2(r))/17 ~ Z' ~ Z. Hence, the subspace 
U characterizes the gauge freedom of the pressure functions. Of course, this can 
also be seen directly. To this end, consider tt S 7^2(0) and q S L 2 {T) such that 
V(7r,9)((^) = 0 for all <p €Y. Then, obviously, tt must be constant in H® for all 
i = 0,... ,k. Integration by parts yields 


(29) 


0 = y [tt] ip ■ i^dx + y 9 Div p dA. 


r r 

For p & Y such that tp • // = 0 on F this reduces to 


0 = y 9 divg p dA, 


proving that q is constant on F® for all i . ,k, and hence [tt] = qH. If F®, 

i = 1,... ,k, is a round sphere, then q is determined on F® by the jump of the bulk 
pressure. On the other hand, if F® is not a round sphere, then |7r] = 0 and 9 = 0. 

Multiplying (1151) with test functions p € Y, integrating over O and F, respec¬ 
tively, and integrating by parts we obtain a weak formulation for this system. For 
functions /i : O —>■ and /a : F —>■ we define 

||/i||-i ,2 := sup fi-pdx and ||/ 3 ||-i ,2 := sup f^-pdA. 


l|V.p|l2<l 


iivv=ii2<i r 


Theorem 3.4. For all (/i,---,/ 4 ) S F^(r) there exists a unique weak solution 
{u,7T,q) G E°(r) of (lisp . We have+++ 

/rJ||V«||2 + /r^||V%||2 + ^||7r||2 + 4ll9l|2 

//2 Li2 

(30) 

< c(—II/ 1 II-I .2 + 4 II/ 2 II 2 + —||/3||_i,2 + 

Mb 

Proof. By Lemma [321 there exists a uq GY such that divMo = (/ 2 , A) and 

(31) ||uo||y < c||(/2,/4)||z- 

Hence, it remains to solve dH]) with /2 = 0, /4 = 0, /i = /i - 2pbdwDuo, and 
fs = fs — ^pDwVuo — [2^b Dut)\v. To this end, we define Fu^ G X' by 

(32) ^“ 0 ( 7 ’)= y {fi-p Y‘2.ph{Duo,Dp)f,) dx + J [fz ■ p F 2p{Vuo,Vp) g) dA, 

n\r r 
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and we let ui G X solve B{ui,ip) = Fuo{(p) for all ip G X. We have 


B{ui,ui) = Fuo{ui) 

< IIVW 1 II 2 II/ 1 II-L 2 + 2fIt\\Dmh\\Xuo\\2 + ||Vui||2||/3||_i,2 

+ 2/r||Pui||2(||VSi;o||2 + ||fc||oo|ko||2) 

(33) w 1 i 1 

< i?(Mi,Ui)2 ( —||/i||-i,2 + (2/16) = IIVuoII2 H-rll-C3ll-i,2 

+ (2/i)^||V®i;o||2 + /i^c_i||Vmo||2), 

where we used the continuity of the trace operator and Poincare’s inequality for 
the second estimate. Combining this with (|3T]) and (ED we obtain 

/Ife l|V“l||2 +M^I|V®Ui|| 2 < c(^||/l||-l,2 + (/l| +C_1 /i5)||/2||2 

■I rll/3ll-i,2 + (ci fj'b + tJ‘^)\\f4\\2)- 

Mb 

Let u := ui -|- uq. Then, from the last inequality and (1311) we obtain the part of 
dSOD estimating the velocities. By Corollary 13.31 there exists a unique (tt, q) G Z 
such that V(7r, q) = B(ui, •) — Fu^- We have 

\B{ui,tp) - Fua{ip)\ 

= J { 2 p-b{Du, Dip)e - fi- (p)dx + j (2/i(T>u, Vip)g - f 3 ■‘p) dA 

n\r r 

< 2/lb||Vu||2||V^||2 + ||/i||- 1,2||V^||2 + ||/3||-|,2l|V(/7||2 

(34) + 2/i(||V®u|| 2 + P|UI|a^||2)(||V®PrMr||2 + IlfclUHMrlb) 

< cr/l=||Vu||2 +/i5||V®u|| 2 + ^||/l||-l,2 + ^ 11 / 311 - 1 , 2 ) 

Mb Mb 

X (m? + M^c_i + -^c-i) ||(/i||v, 

^ Mb" ^ 

where, again, we used the continuity of the trace operator and Poincare’s inequality 
for the second estimate. Note that so far the constants in our estimates do not 
depend on /i//ib. Finally, from the last estimate, the estimate of the velocities, and 
(1281) we obtain (ED- □ 

We saw in the preceding proof that we can easily get estimates with constants 
independent of the Saffman-Delbriick length /i//ib and in which, in particular, we 
can pass to the limit /i 0. Here, however, we are not interested in this limit; cf. 
the remark following the proof of Theorem 13.61 Instead we try to keep the form of 
our estimates as simple as possible. 

In a similar way, we can prove the analogous statement for the system (1181) . 
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Theorem 3.5. For all {fi, - ■ ■, fs) S F;',(r) there exists a unique weak solution 
{u,Tr,q) G IE°(r) of (fT^ . We have 

4l|Vu||2 + M^||V«i;||2 + 4lkll2 + ^lk||2 

(35) < cf"—||/i||_i ,2 + Mh II/ 2 II 2 H—rll/3||-i,2 + M^ll/4||2 

V? h! 

+ -^ll/slb + Mb I/511,2) • 

Proof. We only give a brief sketch of the proof which proceeds analogously to the 
proof of Theorem l3.4l To begin with, let assume that / 2 , / 4 , and /s vanish. Instead 
of X, Y, and Z we consider the spaces 

Xi, := {u G I divu = 0 in n\r, Divu = 0 on T, [u\r G 7J^(r;Tr)}, 

W := {u G | Mr G i?'(r;Tr)}, 

•= {(/2)/4) G L 2 i^) X i 2 (r) such that (fT^ holds with /s = O}. 

Multiplying dT51) by G X^, integrating, and integrating by parts we obtain again 
the weak formulation (1^01) . that is, B{u, (f) = F{ip) for all Lp G X^. Since X^, C X, 
Lemma l3 .1 1 shows that B is coercive on X^, hence giving the unique existence of a 
weak solution u G X^. Analogous to Lemma l3.2l we have to show that the operator 

div : Y„ —>• Z„, (p i-A (div(p, divgMJr) 

is surjective. This, however, is easy to see using L 2 "theory for the Laplacian on T 
and L 2 "theory for the divergence operator in the bulk 17 \r. From this we conclude 
as before that V : —>■ Xf- C Yf, is an isomorphism. In order to deal with full 
data, we construct an extension u' G such that [u']r = /s Then, we let 

u" G Yy solve divu" = (/2 — divu', fi — Divu') and define uq := u' + u". Finally, 
we let iti G Xu solve B{ui,‘p) = Fu^ip) for all p G Xj, for the same as in (1321) . 
reconstruct the pressure functions using the V-operator, and define u := ui+uo- D 

Consider the subspace C L 2 {Pt) x L 2 (r) defined as follows: (tt, g) G Uy iff 
TT = Ki in 17®, TT = kq in g = k® on F® with Kq, k® G M for i = 1,..., fc. Note 
that Zy Z'^'^ {L 2 {^) X L 2 {T))/Uy. Again, the fact that the gauge freedom of the 
pressure functions is characterized by Uy can be seen directly. Indeed, V(7r, q){ip) = 
0 for all Yy implies that n is constant in each connected component of 17 \ F. 
Furthermore, instead of (I29|) . here we obtain 

0 = y g divg if dA, 
r 

proving that g is constant in each connected component of F. 

Now, let us proceed with higher order estimates. Separating the tangential and 
the normal part of by using the equation ([8]) , we see that the system (US can 
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be written in the form 


/ifeAu - gradir = /i 
div u = f 2 

fi[AgV + gTa,dg{w H) + K V — 2diyg{w k)) 

( 36 ) -gradgq + 2fiblDujiy = Prfs 

2fi{{VSv, k)g 2K)) -qH-m = h-v 

divg V — w H = fi 
u = 0 


in ri \ r, 
in ri \ r, 

on r, 
on r, 
on r, 
on do,. 


Note that |Du|u is tangential due to the incompressibility constraint. Indeed, for 
any vector A on F we have [(A • V)u] • u = 0. If A is tangential, we even have 
|(A • V)m] = 0 since u is continuous across F. But then, choosing an orthonormal 
basis ei, 62 at some arbitrary point on F, from div m = 0 we deduce that 

{{v ■ V)u] • V = -|(ei • V)m] • ei - |(e2 • V)u] • 62 = 0. 

Theorem 3.6. Let s > I. For all (/i,..., / 4 ) G F*(F) there exists a unique solution 
(m, 71 ,( 7 ) G ]E®(F) of (fT5)l . We have 

1 ill 

\u\s+l ,2 + k|s-|-l,2 + — kU,2 + —\q\s ,2 

( 37 ) ^ c[—|/l|s-l,2 + Mh l/2k2 + —| 7 ’r/ 3 U-l ,2 + —1/3 • 2 

+ M^|/4k2) +C_s(71fc"||VM||2 +7i^||V«w||2 + + ^Iklb). 

Proof. Theorem 13.41 gives the existence of a unique weak solution (u, tt, q) G E°(F). 
For two reasons we will only give a brief sketch of the proof of regularity. First, the 
general procedure of localization, transformation, and taking difference quotients 
is rather classical. And second, when generalizing the present result to the case of 
a variable membrane, we will have to repeat most of the arguments anyway. 

For integer s we can prove by the techniques just mentioned that our solution 
lies in E®(F) and 

(38) ||(M,T,q)||Eqr) < c ||(/i,..., /4)||F'>(r)- 

By interpolation, the same is true for arbitrary s > I. In order to prove the estimate 
(1321), we only have to replace the lower order parts of norm of the data in (I38|) by the 
terms in the second bracket on the right hand side of (1571) . This, however, follows 
easily by using the equations (1551) and by interpolation and absorption. We present 
the idea by considering the term /i. Let s > 1. For arbitrary integer 0 < 1 < s — 1, 
we have 

I/ll/ < c(MbllV^+‘ull2 + llV^+kk) < c{fVb\u\lU\yu\\l-^o + k|k|k||i-«k- 

Using Young’s inequality and absorbtion, the claim follows. The other terms 
/ 2 ,..., /4 can be handled analogously. 

In the proof of (l38l) we consider only the case s = 1. The case s G N then 
follows by differentiating the equations in the prototype geometries (whole space, 
half space, and double half space). Furthermore, for the sake of a short presentation, 
we assume that (m, tt, g) G E^(F), that is, we only prove a-priori estimates. By a 
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standard localization argument, we may assume that the solution is supported in an 
open cube Qh of side length R> 0. We only consider the case that Qr is centered 
at some point xq G F; the other cases, namely Qr being centered at some point 
cco G on and being contained in 11 \ F, only involve standard analysis of the 
classical Stokes system. Rotating and translating the Cartesian coordinate system 
and choosing R smaller if necessary, we may assume that xq = 0 and that F n Qr 
is the graph of a smooth function h : Q\ := Qr fl x {0}) —>■ (—i?/2, R/2) such 
that h(0) = 0 and V/i(0) = 0. Consider the smooth diffeomorphism 


-Qr^ Qr ■■= '^h^{QR),{x',X'i) {x',X3 - h{x')). 


This diffeomorphism induces the metric e := <i>^e on Qr. We denote the restriction 
of e to Ql by g. Note that : {Qr,Q {QR,e) and ■ {Q%9) (F n 

Qr, g) are isometries. Let us denote the pullbacks of the involved fields by u := 
n := $*7r, v := w := $>, q := fj := ^UPrfs), fs ■= Kih ' and 
fi = for 1 = 1,2,4. By exploiting naturality of covariant differentiation under 
isometries, from (1361) we obtain 


(39) 


Hb^eU - gradg tt = fi 
divg u = f 2 

fi{AgV + gradg(w iFg) + KgV -2 divg(w fcg)) 

- gradg q + 2fj.blD’'ulvs = fj 
2/r((V®u, kg)~g - w {Hi - 2Kg)) - qHg - |if] = 

divg V - w He = fi 

u — V — WVe=0 


in Qr \ Q^, 
in Qr \ Q% 

on Qfij 
on Qfij 
on 

on Q\. 


Here, we included the transformation of the identity u = v + w r on T; the reason 
for this will become clear in a moment. The results from Appendix [B] show that 
(I39|) can be written in the form 


gbAu - gradif = /i 
div u = f 2 

gAv - gradg + 2gblDuje3 = fJ 

=f3 
div V = f 4 
U - V - W 63 = fb 


in Qr \ Q% 
in Qr \ Q\, 
on Q\, 
on Q\, 
on Q\, 
on QIj, 


( 40 ) 
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where 63 := is the normal to x { 0 } and 

A = /i + (e — e) * gradir) -|- r(e) * ((V^e, (Ve)^) * u + Ve * Vm), 

/2 = /2 +r(e) *Ve*u, 

/7 = /7 + (e - e) * gradg q) + ^.b r(e) * ([Vu] + Ve * [u]) 

+ fj, r(e) * ((V^e, (Ve)^) * [m] -|- Ve * [Vm]) , 

A*" = A*" + /rr(e) * (Ve * V®z; -I- (Ve)^ * [rt]) + r(e) * Veg, 
fi = A +7’(e) * Ve * [u], 

/a = (e — e) * r(e) w. 

In fact, the precise form of the functions Ai • ■ • > A given above is not so relevant 
here (except for the terms involving highest order derivatives); it will however be 
crucial when we generalize the present result to the case of a variable membrane. 
It is not hard to see that 


(41) e(a;', cca) — e = r(V/i(x')). 

We may assume without restriction that / 2 , fi, and /s vanish; cf. the proof of 
Lemma 15^ We multiply equation (HiHl i by dfu, where i = 1,2, integrate over Qr, 
and integrate by parts to obtain 


Mb 


\ydiu\ dx + fj, / \'^dii 






Qr\QI 


Qr 


Q 


2 

R 


< llAll2||aff2||2 + |iA^ll2l|aAll2 + c|A^lA9.^li 

< c(||A||2 + l/3-^|i)ll^^*^ll2 + 11 / 3 ^ II 2 II 9 AII 2 . 


For the second inequality we used the Fourier transform and Plancherel’s identity, 
while for the last inequality we used the mapping properties of the trace operator 
and Poincare’s inequality; recall that all involved functions are supported in Qr. 
Using Young’s inequality, absorption, and a standard trick from the theory of the 
classical Stokes system to obtain estimates of all derivatives, see for instance m, 
we infer that 


(42) 


ill 1 

Mb | u | 2.2 + M " | l '| 2,2 + —|| t || 2,2 + —|| g || 2,2 

M? 

< c(^iiAi|2 + ^iiA^ 1I2 + ^1/3^11.2) 

Vb" M^ ^ 


Due to (|4T]l . for R\0, the factor e — e gets arbitrarily small in L°°. Thus, we can 
absorb the highest order terms in A and fj (and A for foU data) on the left hand 
side of (H^ . Transforming back to Qr, summing the finite number of local 
estimates, and eliminating the lower order terms on the right hand side by using 
interpolation and (1501) (or a standard contradiction argument) we obtain (1551) for 
s = 1. □ 


In the preceding proof we took the geometric pullback of the unkowns when 
transforming the system to the prototype geometry. Apart from being very conve¬ 
nient for computations, here, this procedure is indispensible since the normal and 
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the tangential part of the velocity field on T have different regularities. Without 
this geometric transformation we would not be able to control the perturbations in 
fj and /g-L. 

The constants in (1571) depend on ^jl /in particular, the estimate doesn’t hold 
uniformly for ^ 0. The reason for this is not that our proof is too rough; for 

^ = 0 we have to assume Pr/s S (r) in order to obtain u £ \ T). Of 

course, it is possible to establish estimates that hold uniformly for /r \ 0. However, 
since then ([SSIg has to hold in instead of such estimates need 

a higher regularity for the membrane than (1371) : cf. the next subsection. For this 
reason we refrain from giving uniform estimates. 

In a very similar way we can prove the following theorem. 

Theorem 3.7. Let s > 1. For all (/i,..., /g) £ F® (T) there exists a unique solution 
{u,Tr,q) £ lEy(r) of (1181) . We have 



3.2. Analysis for variable membranes. Now, we transfer the results from the 
preceding subsection to the case of variable membranes. More precisely, we will 
analyze the dependence of the constants in the estimates (IdH]) . ([35|), (l37)l . and (|43ll 
on the membrane. We will represent variations of the membrane in the following 
way. Let us denote by S'q,, a > 0, the open set of points in H whose distance from 
r is less than a. It’s a well-known fact from elementary differential geometry that 
there exists a maximal k > 0 such that the mapping 

A : r X (—K, k) —>■ S'k, {x,d) X + d v{x) 

is a diffeomorphism. For functions h : F — >■ (—n, k) we define F?, := {A(a;, h{x)) \ x £ 
F}, and we write x i—>■ {T(x),d(x)) for the inverse mapping A“^. Analogously to 
the notation used before, we denote by F^ and i = l,...,fc, the connected 
components of F/i and the open sets enclosed by F)^, respectively, and by the 
exterior part of \ F/i. 

The analysis will be based on the Hanzawa transform which maps H dif- 
feomorphically to itself such that $/i(r) = F^i. We choose a real-valued function 
/3 £ that is 0 in neighborhoods of —1 and 1, and 1 in a neighborhood of 0 

such that l/3'l < K/\\h\\L^(^r) on F. While, in H \ Sk, we let be the identity, we 
define in 5^ by 

X 1 -^ X + i’(t{x)) h(T{x)) I3{d(x)/ k). 

It is not hard to prove that ^ ^ ^ and ^ph ■= $/t|r : L —F^, are dif- 

feomorphisms; see for instance |22] . Throughout the whole subsection we denote 
by rfi analytic functions of||/i/«:||oo, ||V/i||oo, and, in the case of Theorem 13.131 
and Theorem 13.141 of |h| 2,2 that, additionally, may depend on H, F, p./p.b and s. 
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Furthermore, we write ku, Hh, etc. for the geometric quantities corresponding to 
Th- 

The analysis of the constants in (157)) and (H51) will proceed by a perturbation 
argument. However, this argument inevitably produces lower order terms on the 
right hand sides of the estimates; for this reason the analysis of (|30l) and (l35)) has 
to be done in a different, more direct way0 We will need the following lemmas. 

Lemma 3.8 (Trace operator). Let h G lT^’°°(r) such that |l/i|| 2 ,oo(p) < k/ 2. For 
u G Hq(LI) we have 

IIMr ,,||4 + IMr,,|i,2 < rh\\Vu\\2; 

in particular 

(44) ||MrJ|2<r,|r|i||Vu||2. 

Proof. To begin with, let us assume that h = 0. Then, it is well-know that 

(45) ||Mr||4<c||Vu||2; 
see for instance m- For non-vanishing h we note that 

J\u\‘^dAh = j \uo iph\'^\detdiph\dA<rh j \uo dA 


< 


y |V(uo$^)|2da;^ < rh(^ j \{Vu)o<^jfj\^ dx^ 

Q. O 

iVupl detd($“^)| < rh(^ j\Vda^ . 


n o 

Here, the determinant det d^ph has to be taken with respect to orthonormal bases 
in the respective tangent spaces. For the inequalities we used (1451) and the obvious 
facts that dph and can be bounded in L°° by rh- In particular, we have 

|r?i| < T;,|r|, proving (|44ll . The estimate of the trace in proceeds very 

similarly. □ 

Lemma 3.9 (Divergence equation). Let h G IF^’°°(r) such that ||/i||oo < k/2, and 
let f G ^ 2 ( 14 ) be mean value free in each connected component of Q\Th. Then 
there exists a function u G Ho(H \ F^t) such that divu = / in Ll\Th and 

(46) ||Vu||2<r^||/||2. 

Proof. Theorem HI.3.1 in [16] gives the existence of a solution. This theorem also 
provides an upper bound for the constant in ()46p l^ In order to prove ()46l) we have 
to cover H \ F^i suitably by sets each of which is star-shaped with respect to a 
ball. However, we will only briefly sketch this procedure since the article already is 
quite long. In fact, we solve the divergence equation in each connected component 
of H \ Fft separately. Hence, let us fix some i G {0,..., fc} and consider the set 
We cover Ll\ \ Ss,.. by finitely many balls whose closures are contained in 


®We could also argue by contradiction to prove that ll30t and 113511 hold uniformly as long 
as h remains bounded in a suitable norm. Such an abstract result, however, is insufficient, for 
instance, in the analysis of singularities in the dynamical problems. Furthermore, a Neumann 
series argument would need the assumption of h being small in some high regularity class which, 
again, is insufficient for our purposes. 

®See |13| for an improvement of this bound. 



24 


DANIEL LENGELER 


these balls do not depend on h. Now, let us fix a point x G , let Z denote the 
cylinder of length 2k and radius r that is centered at x and perpendicular to 
and let B C be the ball of radius r/2 centered in Z at a distance 3 k/ 4 from 
r*. We claim that Z is star-shaped with respect to the ball B provided that 
r = l/ci(l + ||Vh|U^ ^(r)) ^ for a sufficiently large constant ci > 0. Indeed, if F* 
is flat in a neighbourhood of x we may choose r = k/ 4(1 -|- ||Vh|| 2 ,^(r))~^, since 
then, as is not hard to see, B is contained in each cone K C with its vertex 
on and its axis perpendicular to T^jF*. On the other hand, if F* is not flat in a 
neighbourhood of x, then we reduce the situation to the flat case with the help of a 
flattening submanifold chart. Even though the cylinders, cones, and balls that we 
construct in this way are not Euclidean, they do contain corresponding Euclidean 
objects with smaller radii and open angles, respectively. This results in multiplying 
r by an additional factor that depends on F* but is independent of h. Now, we 
cover Sljj fl by such cylinders whose inverse radii are bounded by ciTh- Finally, 
in order to control the constants Ck in Lemma III.3.2 in m we have to enumerate 
our covering sets suitably. Here, it is convenient to start the enumeration with the 
cylinders and to end it with the balls that cover \ and that do not depend 
on h] then it is easily possible to control the quantities |Efc| in Lemma III.3.2 from 
below. □ 


Lemma 3.10 (Poincare’s inequality). Let h G 1F^’°°(F) such that ||h||oo < r/2. 
For u G iF^(F/j) we have 

||u - U\\2 < Th Ci||Vu||2, 
where u := |F|j|“^ f u dAh on F), for all i = 1,... ,k. 


Proof. This follows from Cheeger’s inequality, see [IT], stating that the smallest 
positive eigenvalue of the Laplacian on Th is bounded from below by i(Fft,)^/4, 
where 


i{Th) = inf 


ItI 

min(|Hi|,|H2|)' 


Here, the infimum is taken over all simple closed curves 7 dividing F?t into two 
subsets Ai and A 2 with common boundary dAi = dA 2 = 7 , and I 7 I, \Ai\ denote 
the length of 7 and the area of Ai, respectively. From (l69ll we see that is 

bounded by r;, j(F)“^. □ 


We begin with the analysis of the constant in (l30l) . 


Theorem 3.11. Let h G W‘^'°°{T) fl H^{T) such that ||/i||ioo(r) < k/2. For all 
(/i,...,/ 4 ) G F^(F;i) there exists a unique weak solution (u, tt, g) G E°(F;i) of pSp . 
The estimate ioi) holds with the constant c being bounded by 

(47) rh{l + \h\l^oo + \h\l%)(l + Y.i™Hh,Hh)ri - (1,H.)^J“") Cig. 


Here, the sum is taken over all i such that Fj^ is not a round sphere; note that, by 
Jensen’s inequality, the denominators in the second bracket are positive. 


Proof. The proof consists of four steps. We define Xh, Yh, Zh, and Bn just like X, 
Y, Z, and B, respectively, with F replaced by F/i. 
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(a) First, we need to analyze the constant c in (l?T]) 9. From the proof of Lemma IXT] 
and using Lemma [3.81 for u € Xh we obtain 

||Vn||i<^i?,(n,«) + 16||fc,||L||MrJli 

< -Bh{u,u) + rh\\kh\\loCi\\Vu\\l. 

Thus, in view of (l6^ . the definition of Bh, and (1221) . we infer that 
||Vn||^ < -B(w,u)+r,.(l + |/i|^ )ci||Vm||^ 

< rh{l + \h\l^^)—Bhiu,u). 


(b) Next, we analyze the constant c occuring in (|2^ and (|28)) . Let (/ 2 ,/ 4 ) & Zh- 
For each i = 1,..., fc we choose ru = + XW^W'^Hh on such that 

Hh)ri {Hh,Hh)ri = —(1, /4)rj^5 

AUl,l)rl + = (l,/2)ao 

where (•, Orj, (t Onf, denote the L 2 -scalar products on and respectively. 
If F^ is a round sphere, we choose 

AllFj^l :=(l,/ 2 )a^ and A^ := 0. 


If F^ is a not a round sphere, then the linear system is regular, and we obtain 


-ir; 


ll-’ 




■((li/4)r;(l, l)r*, + (lj/2)oj,(l!-ff/t)rj^)- 


Now, let Ip G H'^iTh) be a function with vanishing mean value on solving 
Agi/> = fi+ wHh on F*^, for alH = 1,..., fc. Note that ip.^f = ip-f^ + ip;aK. 
Hence, for v := grad^ ip, integration by parts and (1691) give 


\\XMl = \M\\l- I K^\v\^dA^ 

Th 

< cdlAllz + WwHhWl) +rh{l + |h|2,oo)IMl2- 
However, by Lemma l3.10l we have 

/ |up dAh = - [ifi + wHh)ipdAh < 11/4 + wHh\\2\\ip\\2 
(48) P;. P;. 

< (II/ 4 II 2 + ||w-H/i|| 2 )r,.ci||z;|| 2 . 


Combining the last two inequalities and making use of (1691) we obtain 


||V®u ||2 < Cl rh{l + |h| 2 .oo) (II/ 4 II 2 + \\wHhh) 

<C3r,(l + |h|3_J(||/4||2 + ^|Al| + |A^|). 

i 


( 49 ) 
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Now, we extend the function v + wvh from F/i to Let x : R —^ [0, oo) be a cutoff 
function being 1 in (—k/2,k/ 2) and 0 outside of (—3 k/4, 3k/ 4). Furthermore, let 
ipo be 0 in n \ S'k, and for x € Sk we define 

(po(x) := {v + w iyh){‘Ph{T{x))) xis{x)). 

Then, from (|4^ . (l48l) . and (l6^ we infer that ipo G -ffo(fl) and 

liVtpolb < r/i(K^(||V®n||2 + ||Vw||2 + \\w khh) + i^~^{\\v\\2 + llw^lU)) 

< C7r/i(l + |fr|2,oo + 1 ^ 13 , 2 ) (II/ 4 II 2 + ^ 1-^11 + 1 ^ 2 !) ■ 

i 

Let S iLo(n \ F/i) solve (l27)l in n \ F/j and define := + ^px. From the last 

estimate and Lemma 13.91 we finally obtain 

(50) < r’;i||/ 2||2 + C7r;i(l + \h\l^^ + \h\l 2 ) (II/ 4 II 2 + l^il + l^ 2 l)- 

i 

The fact that this estimate and (|4^ involve the numbers A* is very natural; situa¬ 
tions in which one of the F)^ tends to a sphere without / 2 , /4 fulfilling (asymptoti¬ 
cally) the compatibility condition (I16p have to be penalized. Unfortunately, taking 
the dual of these inequalities doesn’t lead to a suitable estimate for V, at least not 
without further work. For this reason, we simplify the above estimates using the 
fact that 

El^ll + I^2l 

i 

<rh{l + \h\U (1 + E (in|(Lfi„iLi.)r. - ( 1 ,)'') (c 7 II/2II2 + C4II/4II2), 

i 

where the sum on the right hand side is taken over all i such that is not a round 
sphere. Combining this estimate with (|ini) and (1501) we obtain 

11‘^lln 

<rh{l + |fr|E + l^iy (1 + E - (l,i?,.)^J-')cr||(/2,/4)|U.. 

i 

(c) Next, we have to take a closer look at the constants occuring in the proof of 
Theorem (13.41) . Using Lemma [3.81 and (1501) we see that the constant c_i in (ED 
can be bounded by 

rh{l + |fr|2,oo) Cl . 

From this fact and steps (a) and (b) we obtain 

4l|Vni||2 + M^I|V®ni||2 

< ruil + |fr|E + \h\l2) (1 + E - ( 1 , )”') cio 

X (^||/i||-1.2+4||/2||2 + ^11/311-1,2 +M^II/4||2). 

Mb Mb 

Using Lemma Em and (|50)) we see that the constant c_i in (IMl) can be bounded 
by 

T’b(l + |fr| 2 .oo)ci 
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while the constant c_i can be bounded by 

+ I^l2.oo)ci 

and the constant c is independent of h. From these facts, (EU, and steps (a) and 
(b) we obtain 

\B{ui,ip) - Fuoiip)\ 

< rh{l + \h\l% + (l + ^ C 12 



(d) Finally, we combine the above estimates to obtain the stated inequality for 
u = ui — uq and (tt, q) = V“^(i?(ui, •) — F). □ 

As already pointed out in the preceding proof, the occurence of the possibly 
small denominators in (HTl) is related to the fact that the whole expression contains 
no control of the compatibility condition (1161) that / 2 , /4 must fulfill asymptotically 
as one of the tends to a round sphere. If one wants to do analysis in the neigh¬ 
bourhood of a round sphere it might become necessary to prove a refined version of 
m that resolves the singular structure in a more precise way. Furthermore, note 
that (|47l) is not optimal with respect to the regularity requirements on h and with 
respect to the exponents of the semi-norms of h. 

Let us proceed with the analysis of the constant in (l35)l . 

Theorem 3.12. Let h G W^'°°(r) such that ||h||pcx,(p) < k / 2 . For all (/i,..., /i) G 
Fy(r/i) there exists a unique weak solution {u,Tr,q) G E°(r;i) of psp . The estimate 
dSSl) holds with the constant c being bounded by 

rh{l + \h\l^oo) cs- 

Proof. We can proceed very similarly to the proof of Theorem 13.Ill Here, however, 
there is no singular behavior in the vicinity of round spheres since we don’t need 
to prescribe additional compatibility conditions on these. Furthermore, the con¬ 
struction of Uq is much less involved than in the proof of Theorem 13.111 resulting 
in lower regularity requirements on h and in a lower exponent. □ 

Next, we analyze the constants in dSH). 


Theorem 3.13. Let s > 1 and h G nlF^’°°(r) such that ||/i||pcx,(p) < k / 2 . 

For all (/i,...,/ 4 ) G F®(r;i) there exists a unique solution {u,Tr,q) G of 

dnp. The estimate dSH) holds with the constant c being bounded by ru and the 
constant C-s being bounded by 


rh{l + ci\h\l^\h\/l^^2 + Cs+. 


milk 


); 


note that for s > 1 we can get rid of the term involving |h|| by interpolation. 


Proof. It suffices to prove the estimates for smooth h and smooth (it, tt, q). Indeed, 
then, approximating h by smooth functions, applying Theorem 13.61 and Theorem 
13.111 and taking the limit, the claim follows. 

As in the proof of Theorem 13.61 we write u := n := v := 

w := q := ^^q, /* = ^lf^ for i = 1, 2,4, and fs := fj + fs ^3 := ^^(Pr/a) + 
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■ ^) 63 - As before, by naturality of covariant differentiation under isometries, 
we obtain 


- gradgTT =/i inr 2 \r, 
divg u = f 2 in \ r, 
fj.{AgV + gradg{w iJg) + KgV -2 divg(w fcg)) 

(52) - gradgq + 2nblD^^i'e = fj on T, 

2^i{{VH, ke)~g -w{HI- 2Kg)) -qHe-m= ft on T, 

divg V — w He = fi on F, 

U — V — W Ve = fb on r 


with /s = 0 . 


(a) To begin with, let us assume that e is constant in fl and that /s € is non¬ 

vanishing. In this case we can proceed by localization and transformation as in the 
proof of Theorem l3.6l to prove the estimates for integer s and then interpolate for the 
general case. The only difference is that one has to consider a non-Euclidean, but 
constant Riemannian metric in the prototype geometries. In Cartesian coordinates 
the Laplace operators take the form 


Agu = in the bulk, AgV = 5 “^ on Q^. 


In this case, the constants in (l37l) are polynomial functions of the ratio of the largest 
and the smallest eigenvalue of e. 


(b) In the next step, let := e^^(xo) for some fixed xq € fi and assume ||e®^ —e'-f ||oo 
to be sufficiently small. We need to write the right hand sides in (15^ in a more 
precise form. The results from Appendix |B] show that (15^ can be written as 


pibAgu - gradg n = fi 

dive u = f2 

n{AgV + gradg(u) Hg) + KgV -2 divg(zi; kg)) 

- gradg q + 2^ib[D^ulvg = fj 
2^i{{VH, kg)g - w (iJ| - 2Kg)) -qHg-m= 

divg V - w Hg = fi 

u -V -wi/g = f 5 


in n \ r, 

in fl \ r, 

on r, 
on r, 
on r, 
on r 


( 53 ) 
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with 

A = /i + (e - e) * r(e) * gradir) 

-I- fib r{e, e) * ((V^e, (Ve)^) * u -|- Ve * Vu), 

/2 = /2 + r(e, e) * Ve * u, 

fl = /A + (e - e) * r(e) * (/x(V®)^v, gradg q) + fib r(e, e) * {[Vu] + Ve * [u]) 
(54) -I- fir{e, e) * Vfc, fc * Ve, V^e, (Ve)^) * [m] -I- (fc, Ve) * [Vu]), 

/A = /A + r(e, e) * ((fc, Ve) * V®i) + (/c^, /c * Ve, (Ve)^) * [u]) 

-|-r(e, e) * (fc,Ve) g, 

A = A +r(e, e) * (A:, Ve) * [u], 

A = (e - e) * r(e, e) w. 

Here, we replaced V® and V® by V = V® and V® by absorbing the additional 
coefficients into r(e) and r(e, e), respectively. It is not hard to see that in 5”^ 

e — e = r(h/K, hk, Vh). 

Here and in the following lines, for simplicity of notation, we write h and k instead 
of ft.or and k or. From this identity, by induction we infer that for I G N in Sk we 
have 


l-\-l 


= i^^^r)ih/K,hk,vh)hP'^ *Y[iiv^yh)p* *Y[i{v^ykY 


.-p 


P,Pi^Pi>0 




where stands for derivatives of r with respect to the spatial variable and 

/+i i i 

-Po+Pi “ 1)k + + 1)A + P = l, Po < + 1)A- 

z— 2 z— 0 z— 0 

In particular, we have ~ ^)Pi — From this we easily deduce that in 


(55) 


i 

VV(e, e) = r(e, e) * t)[(V^e)^^ 

1j>0 j=i 

i+i 

= ^r(ft/R,ftft,Vft,e)*n((^A*/rr, 

Pi>0 i—2 


where Xj=i ^ ^rid XAA* ~ ^)Pi — Now, fix some S G [0,1). For ft > 0 and 
a G {I + ft, 2, 2 +ft,...} let 


la {1 + ft, 2, 2 -|- ft,..., a}. 


while for ft = 0 and a G N >2 let 

la ■■= {2,3,...,a}. 

Finally, for a as above and 6, c G {0, ft, 1,1 + ft,...} let 

Ja‘' '■= [iPiJ)ieia M G {0,ft, I,l + ft, ...,c} and 

Pi G N>o such that ^ {i — l)pi + I <b 
i&Ia 
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(b.l) We start with the analysis of the term which is, with regard to the expo¬ 
nent, the limiting one. We write s-|-^=n-|-(5fornG N>i and 6 G [0,1). Starting 
from the expression for in (|M)) and using (1551) , Lemma 1C.21 Corollary 1C.31 and 
we infer that 


~rl/3'ls-i,2 ^ ~i7^l/3"L-i.2 
Mb Mb 


ieM+ 3/2 


Mb 


(Pi 


1 \C 1 7» + 3 / 2 .s- 1/2 
3 + 3/2 


Th 


n I'-i 


i'piri,2\m\l,r,2 


lei, 


s + 3/2 


+ ^ ’''*11 \’^\^Pir.,2\Q\l,r,2, 

»6L+3/2 


where 2 < r, < oo are chosen such that 


1 

r 


+ E 

i-. Pi^O 


1 


PiTi 


1 

2' 


Here and in the following, with a slight abuse of notation we define | • \o,r ,2 ■= II' ||r 
for ease of notation. 

To begin with, let us consider the term involving |'5|s+i/2,r,2 and p 2 = 1; in this 
case all other pi vanish since — l)pi < 1. Assuming 2 < r < 4, we use (1711) 
with t = 2, p = r 2 = 2r/(r — 2), and a = 2,/{sr) as well as (1731) with t = s + 1/2, 
p = r, and /3 = (s -|- 1 — 2/r)/(s -|- 1/2); for the first estimate we have to exclude 
the case r = 2 if s = 1. Hence, we have 


-jWl^|2,r2.2|fl|s+i r.2 < 

Mb 


Th 





The exponent 

— _= ^2 -I- _ - _ 

1-/3 V s/2-r/2 

of |/i|s_|_ 2,2 is minimized for r = 2. Furthermore, e has to be chosen so small that 
the |'i;|s+i_ 2 -term can be absorbed on the left hand side of the estimate proved in 
step (a). However, the constants of continuity in the latter estimate are of the form 
Th; hence, e can be chosen to have the form Thus, for s > 1 we have 


(56) 


M 

1/2 

Mb 


\h\ 


2,oo|'y|s+i 2 ^ 


< rh(l + (^) l^lE/2)Mb l|Vu||2 + e/x3 


^^|s-|-1.2- 


For s = 1 we use the interpolation estimate only for v. This gives 


M 

1/2 I 
Mb 


|/j|2,oo|fi||,2 < + (^) |/*l|oo)Mb ll^^lb + eM" I'5|3-11.2. 


Next, let us consider the remaining terms involving |il|s+i/2,r,2- Note that in this 
case all pi with i > 2 vanish. Again, we choose r = 2 and = oo for all i and 
consider (jTTl) with t = i, p = oo, and a = 9i := {i — l)/s; note that 9i G (0,1). 
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From this estimate and (17^ with t = s + 1/2, p = 2, and j3 = s/{s-\-l /2) we obtain 


1/2 

^^b ieIs+3/2 


n i^Iloo, 




2 < 


Since ^iPi — ^ this gives 


n i^Iloo.2 


Pb 

+ efj,^\v\s+i, 2 - 

, S+l 


It'L+i 2 < 


rhU + (—) l|Vw ||2 + e/^" 


iih+1.2- 


1/2 

Pb ieis + 3/2 

Now, we assume that d > 0 and consider the terms involving |u|„^r, 2 - In this 
case we have ^(/ — l)pi < 1 + d and hence all pi with i > 3 vanish. For s > 3/2, 
we can choose r = 2 and Vi = oo for all i since then we have Oi S (0,1) even 
for / = 2 + 5. Using GD) with t = i and a = as well as (iza) with t = n and 
/3 = (n — l/2)/(s -I- 1/2) we obtain 


(57) 


1/2 

‘i‘^Is+3/2 


n l^li'./o.2|i^l».2 < r,*(l+ ( —) "|Il|E2)2’)^/l6 llVulh 


Kflb' 

+ e^5|i;|s+i,2- 


Since ^iPi ^ (1 + ^nd 
1 + S 1 


1-/9 


(-0 


1\ 1 + 5 


sJ 1 + 26 


<2 


we can dominate the right hand side of (l57l) by the right hand side of (1561) . However, 
in the case s < 3/2 we need to proceed differently. In this case we have n = 1 and 
<5 G [1/2,1). We choose r = 4, rj = 4, and = oo for / j; here and in the 
following we denote by j the largest i G 7s-i-3/2 such that pi ^ 0. We want to use 
(1711) with t = j, p = Apj, and a = 0' := (/ — 1 — l/(2pj))/s. Note that 

d 


0 < < ,; < FG 


= 1 ; 


S ■' S 

in the case <5 = 1/2 we may assume without restriction that pj > 1 so that 0' >0. 
Using this estimate as well as (GU with t = i, p = oo, and a = Oi and G1 with 
t = n = 1, p = A, and P = l/{s + 1/2) we obtain 

n \^\^!oo,2\h\^^Ap,,2\v\lA 

(58) 

<’'/*(!+ (^) " '”V?l|Vu ||2 + ep^\v\s+l, 2 - 

Since ^ OiPi + O'.pj < (1/2 + 5)/s and 


1/2 + 6 1 




s 1 — /3 

for s > 1 we can dominate the right hand side of (1551) by the right hand side of 
(1551) . 

Now, we consider the terms involving |[u]|s_i ^ 2 - In this case we have X)(* — 

l)Pi < 2; however, since //“ only contains up to first derivatives (squared) of e 
and no second derivatives, in fact all pi with i > 2 vanish. We choose r = 2 and 
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Vi = oo for all i and use jni) with t = i, p = oo, and a = 9i as well as (ESI) with 
t = s — 1/2, p = 2 , and 7 = (s — l)/s; again, we have to exclude the case s = 1 so 
that 62 = l/s G (0,1). This gives 



1 

+ ep^\u\s+i, 2 - 


Since ^ 6 iPi < 2/s and 1/(1 — 7 ) = s we can infer that 

|/iir;oo.2lNls-i,2.2 < (—) l^ls+2,2)^6 l|Vu||2 

(59) ^ ^ 

1 

+ |u|s+i.2- 

For s = 1 we simply don’t use the estimate (EID for t = 2 and obtain 


K . = V ^ 

1 

+ eMb 

An application of Young’s inequality yields 

-^rh Jl l^lLoo,2lMli,2,2 < + —(I^l3,2 + I*I2 ,oo))KI|Vu||2 

( 60 ) ^6 * 6 A /2 

1 

+ £^6 |w| 2 , 2 - 


Next, we assume that <5 > 0 and consider the terms involving |['S]|n-i,r-, 2 - In this 
case we have ~ ^)Pi ^ 2 + (5; however, all pi with z > 3 vanish. As before, we 
have to distinguish the cases s > 3/2 and s < 3/2. In the case s > 3/2 we choose 
r = 2 and Vi = 00 ; note that 9i S (0,1) even for i = 2 + 5. We apply (Ell) with t = i 
and a = 9i and ESD with t = n—1 and 7 = (n — 3/2)/s and proceed as in the case 
I = s — 1/2. Since — (2 + i5)/s and 

2 + ^ I _ 2 + 5 ^ 
s 1 — /3 1 + 5“ 


we can dominate our terms by the right hand sides of (15^ and (IBUl) . respectively. 
In the case s < 3/2 we note that 5 > 1/2 and choose r = 4, rj = 4, and = 00 for 
all i ^ j. Note that i < 2 and hence 9i G (0,1), except for 6*2 in the case s = I; this 
case has to be treated separately, exactly as above. Furthermore, note that 

1/2 + 5 


^-l/(2p,-) < < 1/2 + 


< 1 ; 


S ■' s 

without loss we may assume pj to be so large that 9/ > 0. Now, we apply (1711) 
with t = i and a = 9i, (1711) with t = j and a = 9/, and (1751) with t = n — \ and 
7 = (n—I)/s and proceed as in the case I = s—1/2. Since (^iPi+9/pj < (3/2+5)/s 
and 


3/2 + 5 1 3/2 + 5 ^ ^ 

s I-/3 “ 1/2 + 5 “ 

we can dominate our terms by the right hand sides of (l59l) and (I60p . respectively. 

Next, we consider the terms involving | [it] |i,r -,2 with 0 < I < s — 3/2. In this 
case we have ^ s + 3/2 — /; however, all pi with i > s + 3/2 vanish. 
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We choose r = 4, rj = 4, and Vi = oo for all i ^ j. Note that for all i ^ j we 
have i < n + 1 in the case 6 > 0 and i < n in the case <5 = 0; therefore, we have 
9i G (0,1). Furthermore, we note that Ps+ 3/2 < 1 so that 

of ^ ^ ~ ^ 1 . 

^s+ 3/2 < -^< 1 . 

similarly we see that 9!^^^ < 1, and 0' < 1 holds trivially for j < s + 1/2. We apply 
(ITT]) with t = i and a = 9i, (1711) with t = j and a = 9j, and (1751) with t = I and 
') = l/s. Since ^ 9iPi -\- 9jPj < {s — I + l)/s and 

s — ^“t“l 1 s — / “t“l 5 

s 1 — /3 s — I “3 

we can dominate our terms by the right hand sides of (15911 and (I60L respectively. 

The terms involving the surface pressure q can be treated almost exactly like 
the terms containing the trace [u] of the bulk velocity; the main difference is that 
one has to interpolate between |g|s _2 and II 9 II 2 instead of |m|s+i ,2 and ||Vm|| 2 . We 
only present the case I = s — 1/2 where — l)pi < 1. We choose r = 2 and 
ri = 00 for all z; again we have to treat the case s = 1 separately. Applying (1721) 
with t = s — 1/2 and (I7ip with t = i and a = 9i we obtain 

n \f^\^!oc,2\9\s-l/2,2 <rhU+( — ) \h\'l+2,2)^Uh 

1 

+ e—|u| 2 , 2 - 

/r2 


(b.2) Next, we analyze the term /s. Let s+^=n + (5 for n£ N>i and 6 G [0,1). 
We have 


Mh l/sl 


S+: 


,2 — 


looThP/\U\s+l ,2 


P/Th 


E 

s + l/2,s 

(pi,0^‘^s+3/2 


n 


\hZ.r .,2 


|[fi]|^,r,2 “b I ^1 l-t-(5,ri+5 ,2 I [fi] |n,r,2 


iel. 


+ 3/2 


where b{ 6 ) = 1 if d > 0 and b{ 6 ) = 0 if d = 0 , and where 2 < r, r/ < 00 are chosen 
as above. 

First, let us assume that i5 > 0 and consider the terms involving In this 

case we have ~ ^)Pi — '^i hence, all pi with z > 2 vanish. Choosing r = 2 and 
ri = 00 for all z and applying GH) with t = i and a = 9i and dZSl) with t = n and 
7 = (rz — 1 / 2 )/s we obtain 


(61) 


p/rh |h|f;p,^,_2l[h]|n.2 < T,i(l + p/\\Vu\\2 + ep/\u\s+i:. 

ieIa + 3/2 


<rh{l + \h\s+ 2 , 2 )p/ ||Vzz ||2 + ep/ |uU+i. 2 , 


since Y/ ^ ^/s and 

(5 1 6 s 

si — 7 ss — n+ 1/2 


Next, we consider the terms involving |[zi]|/,T -,2 with 0 < I < s — 1/2. In this 
case we have ~ ^)Pi < s + 1/2 — I- We choose r = 4, rj = 4, and r/ = 


00 
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for all i j and apply (EU with t = i and a = 9i, (EU with t = j and a = 9'^, 
and d75l) with t = I and 7 = l/s. We see exactly like in the analysis of the 
terms involving |[{t]Er .2 with 0 < I < s — 3/2 that 9i G (0,1) and 0' < 1. Since 
Y. (^iPi + (^jPj < (s - 0 /s and 

s — I 1 ^ 


s 1 — /3 

we can dominate our terms by the right hand side of (I 6 ip . 

(b.3) Now, we consider the terms fj and / 4 . Let s = n + ^ with S G [0,1). We 
have 




rl/ 3 '^ |s-1.2 + \fi\s,2 < —IfJ |s-l,2 + \h\s,2 

2 ^2 


+ \\e - e\\oorh{p^^ l'yU+1,2 + ^|g|s, 2 ) 

P2 

+ p'^Th E n l^ltpin,2l[“]kr,2 +0(^)|h|l+5,ri+,,2|ii|„+l.r,2 


+ —Th 

H2 

Pb 




E n l^lLPiG,2l9kD2 + feWI^|l+5.ri+,,2|g|n,r.2^ 


+ 4 E 


fh n l^irPir..2lMkD2, 

(Pt,0GJr+l *6L+l 

where b{5) and 2 < r, < 00 are chosen as above. 

We begin by assuming 5 > 0 and analyze the term involving |'0|ra+i,r,2. Choosing 
r = 2 and ri+s = 00 and applying EU with t — 1 “ 1 “ (5 cind Ck — ^i+(5 Bjiid (jT3[) with 
t = n + 1 and /3 = (s + 1/2 — ( 5)/(5 + 1/2) we obtain 


since 


1 { { 2 1 —|— ^ \ . 1 . 1 

^5r,i|/i|i+5,oo,2kk+i,2 <rh{l+ y—j \h\s+i" 2 )Pb l|Vw ||2 + efj, 0k.+i.2, 

9i+s 


1-/3 ^”^25' 

Next, we consider the terms involving /r 5 |[u]E^ 2 - In this case we have — 
l)Pi < 1 and hence all pi with i > 2 vanish. We choose r = 2 and = oo for all i; 
again, the case s = 1 and P 2 = 1 has to be treated separately. For s > 1 we apply 
d7T|) with t = i and a = 9i and (1751) with t = s and 7 = (s — \l2)ls and obtain 


, 1/2 


n i^iE. 


2INU.2.2 


ie/s 


(62) 

VI 

^1 + 1 

[-) 



< 



VI 

/ +1 

(-) 






ls+ 1.2 


|s+l,2- 


Here, the last inequality follows from Y, (^iPi < 1/s and 

1 1 


s 1 — 7 


= 2 . 
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We omit the case s = 1. 

Now, we assume <5 > 0 and consider the terms involving /x 2 In this case 

we have < l + J and hence all pi with i > 3 vanish. For s > 2 or p 2 +s = 0 

we can again choose r = 2 and = oo for all i since then we have di £ (0,1). We 
use (ITTIl with t = i and a = 6i and (l75|l with t = n and 7 = (n — ll2)ls. Since 
< (1 + <5)/s and 

s l-P 1/2 + 5 - 

we can dominate our terms by the right hand side of (1621) . For s < 2 and p 2 +s = 1 
we choose r = min(2/(1 — 5), 3) and r 2 +s = 2?"/(r — 2). Using (jTTl) with t = 2 + 5, 
p = r' 2 + 5 , and a = (5 + 2/r)/s and (1751) with t = n, p = r, and 7 = (n-|-l /2 — 2 /r)/s, 
we can dominate our term by the right hand side of (1621) since 

Q ^ S + 2/r ^ 2 

1-7 5 + 2 /r-l/ 2 - 

Next, we consider the terms involving |['S]|z,r ,2 with 0 < I < n. In this case 
we choose r = 00 , rj = 2, and Vi = 00 for all i ^ j. We use dni with t = i and 
a = 9i, (1711) with t = j and a = 9” := {j — 1 — l/pj)/s, and (1751) with t = I 

and 7 = (1 + l/2)/s. Note that 7 G (0,1) since I < s — 1, and 9'/ < 1 since 

Ps+ 2 ,Pn +2 G {0,1}. Since X) + (^'/Pj < (s - 0/s and 

s — I 1 s — I ^ 

s 1 — 7 s — I — 1/2 — 

we can dominate our terms by the right hand side of (|62l) . 

The terms involving q can be handled analogously; we only present the case 
I = n. Using (17T1) with t = i, p = 00 , and a = 9i and (1721) with t = n and /3 = n/s 
we obtain 

—r^/i|^|l-|-<5,oo,2|9|n,2 < T/j(l + | h| s_|_ 27 ) —^ || g|| 2 + ^^t\q\s,2 

^2 p ,2 p ,2 

since 

9i+s _ 5 s _ ^ 

1 — P s s — n 

Next, we consider the terms involving p-b/P^[u]\i,r, 2 - We choose r = 3, rj = 6 , 
and Ti = 00 for all i ^ j and apply CD) with t = i and a = 9i, CD with t = j and 
a = 9'J' := (j — I — l/{3pj))/s, and (1751) with t = I and ') = {I — I/6)/s. This gives 

n \Kip,r,,2M\l,r,2<rh{l + 

1 

+ epl |uO-Hi,2. 

Since ^ 9iPi + 9'J'pj < {s — I — l/3)/s and 

s-l-1/3 1 _s-^-I/3^^ 

s 1 — 7 s — ^ + 1 / 6 “ 

we can dominate our terms by the right hand side of (1571) . 


m 




)‘'7jl|vi||! 
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(b.4) It remains to analyse the terms fi and / 2 . Let s = n + S with n € N>i and 
d G [0,1). We have 


— |/l|s-l,2 + fJ-h \f2\s,2 < —|/l|s-l,2 + Mh \f2\s,2 

Mb Mb 

+ ||e - ellooi'bftib l^ls+i.2 H-rl^ls. 2 ) 

^ M? ^ 

+ M?’’/if n '^^~\^f^!piri,2\^\l,r,2 + KS)K'^^+^\h\l+S,r^+s,2\u\n+l,r 


+ Y['^’''\f^f^!pir.,2\^\l.r,2 + HS)^i"^+^\h\l+S,r,+i,2\^\n,r\ 


where b{S) and 2 < r,ri < 00 are chosen as above. Here, we can proceed analogously 
to (b.3). The main difference is that we have to apply the interpolation estimate 
(Tfil) for a = Vti, TT. 

To begin with, let us consider the terms involving |u|i,r ,2 with / = s if <5 = 0 and 
I G {n, s, n + 1} if d > 0. We choose r = 2 and = 00 for all i except for the case 
I = n,s < 2,P2 +s = 1 which has to treated separately as above. Excluding this 
case we apply (fTTll with t = i and a = 9i and (1741) with t = I — 1 and 7 = (^ — l)/s 
and obtain 


Mb^b 


n 




'l^ltp.n.2|wkD2 < rb(l + K"|M|^2*2*) ' "'Mb I|Vm||i 


+ eMb kk+1.2- 


Since ^iPi ^ (s + 1 — O/'S and 

s — Z + 1 1 ^ 

s 1 — 7 

we can dominate our terms by the right hand side of (1621) . In the special case we 
apply (1741) with t = n —1 = 0 , p = r:= min( 2/(1 — d), 3), and 7 = (3/2 — 3/r)/s 
and dZlD with t = 2 + 5, p = 2r/{r — 2), and a= {6 + 2/r)/s. Since 

« ^ 5 + 2lr ^ 2 

1 —7 (5 + 3/r—1/2“ 

we can again dominate our term by the right hand side of (1621) . 

Next, we consider the terms involving \u\i^r ,2 with 0 < Z < s — 1. We choose 
r = 00, Tj = 2 , and = 00 for all i ^ j and apply (1711) with t = i and a = 0i, 
d7T|) with t = j and a = 9"', and (1741) with t = I — \ and 7 = (Z + l/ 2 )/s. Since 
Y. (^iPi + (^'j'Pj < {s-l)/s and 

s — I 1 s — I ^ ^ 

s l-P ~ s-l-l/2 - 

we can again dominate our term by the right hand side of (|62l) . 

Estimating the terms involving the bulk pressure is proceeds completely analo¬ 
gously. 

We want to absorb the highest order terms on the right hand sides with a factor 
||e — e||oo on the left hand side of the estimate proven in step (a). The constant in 
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the latter estimate is of the form rn, so we need to assume ||e — e||oo to be of the 
form . In this case the lemma is proved. 

(c) Now, let us consider the general case. For xq G and i? > 0 we have 

||e - e(a;o)||oo,BH(a:o) < -R||Ve||oo.n < Rrh{l + ||V^/i||oo). 

Since we want to apply step (b), we choose R to be of the form 1| V^/i|loo))~^ 

and cover fi by balls of radius i?/ 2 ; the number of balls we need to do so is of the 
order 

^ ^ = diam(fl) rh (1 + || V^/i||oo). 

Let be a smooth function that is 1 in i?i(0) and 0 in Now, we fix one of the 

finitely many balls B}i(xo) and multiply our solution with ipn := — xq)/R). The 

product is again a solution, however, with right hand sides perturbed by lower order 
derivatives of the original solution multiplied by suitable powers of 1 /R. These lower 
order terms can again be handled by interpolation and absorption. For dimensional 
reasons this must lead to an additional factor of the form t?i (1 + || V^h||^) in front 
of the lower order terms in our estimate. 

(d) So far, we proved the estimate 



We have to transform this estimate back to the original domain. Equivalently, we 
can replace the Euclidean norms by norms which are defined with respect to the 
perturbed metric e; we denote such norms with an index e. We begin by replacing 
the norms of the data. Erom (l66l) and (l76l) as well as Lemma 1C.21 and Corollary 
1C.31 we infer that 


\fl\s-l, 2 <rh( n ^"H^l?p.n. 2 l/lkc 2 

+ 6((5)k’’i+' 5 |/l|i+yri+j,2|/l|rt-l,r + e|/l|s-l,2 ) , 


where s = n + 5 with n G N>i, 6 G [0,1), and b{6) is defined as above. Now, in the 
terms involving \fi\i,r with Z < s — 1 we replace /i by the left hand side of (I52l) i and 
take into account the representation of this left hand side derived in (I551) i (with e 
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replaced by e). This way we see that 

|/l|s-l,2 < fh g|/l|s-l,2 

+ ^J'brh( n 2 +b{6)K'-^+>>\h\i+S,ri+s, 2^71+1,T 

+ —'^h( ,21^^11,r,2+b{6)K'^\h\i+s,r,+sM^\n,r]- 


These are exactly the terms we dealt with in (b.4). While in (b.4) the highest order 
derivatives were absorbed on the left hand side, here, these terms are implicitly con¬ 
tained in the expression e\fi\s-i,r- The terms involving / 2 , • ■ •, A can be handled 
completely analogously. We obtain 


(63) 


1 ill 

|■u|s+1.2 + l^'U+1.2 + —|7rU,2 + —\q\s,2 

< r/i(^e|/l|s-1.2 +M?e'|/2k2 + \s\fj U-1.2 + ^e~|/kL-l 2 

+ e\f4\s,2^ + rii(l + k|2,oo + ^s+f kls-l- 2 , 2 ) 

X (Ai|l|Vu||2 +l^^||V®{i||2 + ^|k||2 + ^Ikk)- 


Next, we take care of the terms on the left hand side. Using again (IMl) and GlD 
we prove that 


|^^U+l .2 < rh 


n klkpin. 2 kkD 2 - 




Again, these are exactly the terms we dealt with in (b.4); here, however the highest 
order derivative is estimated using (1551) instead of being absorbed. The term in¬ 
volving TT can be handled analogously using (b.4), while the terms involving v and 
q can be dealt with using (b.3). From these arguments and by a transformation 
of the left hand side and the data terms to the original domain 11 \ T^ and T/j, 
respectively, we obtain 



It remains to treat the lower order terms on the right hand side. Obviously, for the 
pressure terms we have 



ON A STOKES-TYPE SYSTEM ARISING IN FLUID VESICLE DYNAMICS 


39 


Concerning the velocities, we compute using (1661) 

||Vu ||2 < Th g||V®'u ||2 +rh{l + K^|/i| 2 , 3 )||w ||6 

< r/ig||VS||2 +r/i(l + K5|/l|i^2,2)l|V“l|2, 

where we used dD and Sobolev’s embedding in for the second inequality, and 

/i^||V®u||2 < fJ-^rhg\\V^v\\2 +rh{l + \h\2,i)^Ji^\\v\\i 

1 

< ^J-^rhg\\S/^v\\2+rh[l + ^|/i||^2.2)m?I|Vu||2, 

where we used (ED and Lemma 13.81 for the second inequality. Here, we avoided 
absorption since this would raise the exponent of |h|s_|_ 2_2 unnecessarily. 

This, finally, finishes the proof. □ 

At first sight, it might seem easier to prove the preceding theorem for integer s 
first and apply interpolation afterwards. However, first, this approach would slighty 
worsen the bound on the constant c_s. And second, this would not avoid the need 
for estimating products of functions in | • \t ^2 for non-integer t since the equations 
(I52D4.6 must be analyzed in iJ® ^T) and iJ®+^r), respectively. 

Finally, we analyze the constants in (l43ll . 

Theorem 3.14. Let s > 1 and h G such that ||/i||^oo(p) < k/ 2. For all 

(/i, ..., fi) € FKFfi) there exists a unique solution (u,7r, g) S E®(r/j) of dT51) . The 
estimate (03) holds with the constant c being bounded by rh and the constant C-s 
being bounded by 

r/.(l + Cl |/l|^_^|/l||;2.2 + Cs-pi|/j|s+|;2); 

note that for s > 1 the we can get rid of the term involving |/i|| ^ by interpolation. 

Proof. We can proceed exactly like in the proof of Theorem l3.13l Here, however, we 
don’t need step (b.l), leading to a slightly lower exponent, and step (b. 2 ) becomes 
a bit simpler since w is prescribed and, in particular, absorption is not needed. □ 


Appendix A. Mechanics of fluid interfaces 


In this appendix, for the sake of completeness, we repeat the basic facts con¬ 
cerning the mechanics of fluid interfaces; more detailed information can be found 
for instance in [MlIIlITO]. 

We begin with the general kinematics of fluid interfaces. Let Ft C be an 
interface moving with velocity u = v+w u, where v is tangential and u is a fixed unit 
normal field on Ft, and let g = gt he the Riemannian metric on F( induced by the 
ambient Euclidean space. Consider the familiy of diffeomorphisms tpt.s : Lt —^ Fg 
associated with the vector field u, that is, g)t,t = idpj and dgipt.s = u o tpt.s- We 
denote by Df /Dt the material derivative of a scalar field / on F( with respect to 
the vector held it, that is, 


£1 

Dt 


■— ds\s—tf O 


Throughout this appendix we shall work in convected coordinates (a;“) on F*, that 
is, x°'\t = x“|s o let (da) be the associated coordinate vector helds. Note that 

.-1 /o \ „ ,.-i 


^a|s — {^d(pt^sb)a\t) O Tt.s — O Tt.s 
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and thus, by Schwarz’ theorem, 

~ ^s\s—tda‘^t,s — daU\ 

here, we take the material derivative of each Cartesian component of the vector 
field da- Analogously to the bulk case, the surface rate-of-strain tensor, T)u, is given 
by the rate of change of (infinitesimal) lengths along the flow or, in other words, by 
the Lie derivative of the Riemannian metric. More precisely, if the tangent vector 
field Xt on Ft is transported by the flow, that is, Xs = {dipt,sXt) o then 

Vut{Xt,Xt) ■.= l^^g^{Xt,Xt), 

where | • | denotes the Euclidean norm. In convected coordinates, we have 

mi - 1 ^ D 

{Vu)^0 - 2 “ 2 Dt 

= ^{da{v -\- w n) • 9/3 + • djs{v-\-wn)) 

= ^ (l’/3;a + W l^,a • 9/3 + Va-jS + W I^,/3'9a) 

= ^(v/3;a + Va-,is) -wka/S, 

where we used the identity ka/s = —v,a ■ dp for the third equality. Note that the 
material derivative of the components in convected coordinates of any tangential 
tensor field yields the Lie derivative of this tensor field along the flow; however, 
this is not true for hybrid tensor fields which are partly expressed in Cartesian 
components, see below. Since g°‘^gp^y = is a constant we have 

Again, since = 1 and v ■ da = Q are constants we have 

V^V = 0 and da ■ = -V ■ ^da = -{v ■ daV + w^a) = -kapv^ - w^a- 

From the chain rule daf\s o ^^ 1,3 = da{f o g^t,s) we see that Dt and da commute, 
and thus 

9/3 • = ^/3 ■ ^ -{ka-fV'^ + W,a);/3- 

Hence, we have 


^kad — ^ 


-rrdaV — OaV ■ 

Dt Dt 


= W-aP + {ka-fV'^)-p + kfv^.p - wka-fkj. 


From this we see that the material derivative of the mean curvature H = g°^^kap 
reads 

= kap^g'^^ + g°‘'^^kap = Agw + u>|fcp + k^p^l^ 

= AgW + w{H^ — 2K) + dH{v). 

Here, we used the fact that the tensor kap--y is symmetric with respect to all indices, 
which is the content of the Codazzi-Mainardi equations. Finally, we compute the 
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material derivative of the volume element dA = dAt on F*. To this end, we note 
that 


‘Pt,l{dA,) = 


Vdet Qs o if t,^ 


dAt 


\/det gt 

where the determinants are taken with respect to convected coordinates. Hence, 
we compute 

9s|s=t\/det gs o '^t,s = V'det gt = Vdet gt (gt)“^(^Wt)a/3 


and obtain 


— dA = (divg V — w H) dA = Div u dA. 


Next, let us have a look at conservation of mass and linear momentum on a fluid 
interface. Consider a mass density p on Ft that is being conserved and transported 
by the vector field u; this means that for all subsets Mt C Fj moving along the flow 
we have 

4 / pdA = Q. 


dt 


Mt 


Interchanging differentiation in time and integration in space gives 

Dp 


Ttj + 

Mt Mt 

since Mt is arbitrary, this gives conservation of mass in differential form 

Dp 


Dt 


+ p Div M = 0 on Ft. 


Again, let Mt C Ff be (a smooth bounded domain) moving along the flow. Then, 
for i = 1,2,3, balance of linear momentum in integral form reads 

^ puMA= [ T^ii^9Mtrds + [ fdA, 


(64) 


dt 


Mt 


dMt 


Mt 


where T is the surface stress tensor, vgMt is the outer unit normal to Mt, and / is 
a given surface force density. Gauss’ theorem gives 

I Tf{v9Mtrds = I g^^Tl^pdA. 

dMt Mt 


The integrand on the right hand side, however, is nothing but the f-th component 
of the surface divergence Div T. Interchanging differentiation and integration on 
the left hand side of ((Ml) gives 

S / P^.dA= I (^u + p^ + puDiyu)dA= I p^dA. 

Mt Mt Mt 


where we used conservation of mass for the second identity. Since Mt is arbitrary, 
we obtain conservation of linear momentum in differential form 

= DivT + Z on Fj. 


It is interesting to note that this expression has to be interpreted with respect to 
Cartesian coordinates in since, as we already pointed out, the material derivative 
must be understood component-wise. In other words, while Du/Dt is the Lie 
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derivative of each Cartesian component it* along the flow, it is not the Lie derivative 
of the vector field u (which would vanish anyway due to skew-symmetry of the Lie 
bracket). In particular, balance of linear momentum is not a tensorial postulate; 

cf. [25] . 

Finally, let us compute the stresses resulting from the Canham-Helfrich energy 
O; this can be done by exploiting translational symmetry of the energy. Consider 
a fluid interface Ft that is being translated, that is, (pt,six) = x + ait — s) for some 
a € Then, we have = g°^^a ■ dp and w = a ■ u. Since the energy (|T|) is 
invariant under translations, for all translated subsets Mt C Ft we have 

Mt 

= At / (^^iH — Co)i^AgW + wiH^ — 2K) + dHiv)) 

Mt 

-f (divg v-wH)'^dA. 

Note that the integrand can be written in the form 
{AgH + HiH^/2 - 2K) + Coi2K - HCo/2))w 
+ divg iiH - Co) gradg w -wgradgiH - Co) -b u(7L - Co)^/2) 

= ^ o*(K“^gradi 2 FC -b divg ((iJ - Co)gradg C - Cgradg(iL - Co) 
i 

+ (iL-Co)V2 5“^9^a„). 

Since Mt and a are arbitrary, we conclude that the L 2 -gradient of the Canham- 
Helfrich energy can be written as the surface divergence of the stress tensor 

= At( - (H - Co)fcf5^ -iH- Co).oC + (H - Co)V2 5;), 

that is 

(65) — grad ^2 F v = Div ^T. 

While this formula looks very natural from a physical point of view, from a math¬ 
ematical point of view it is not obvious at all that the L 2 -gradient can be written 
in divergence form. In fact, when this and related identities were discovered in |33] 
by purely mathematical (and rather complicated) considerations, it was a major 
advance in the analysis of Willmore surfaces. Finally, let us verify the identity (l65l) 
directly: 

K-i Div = - {iH - Co)k°‘^) .^dp - iH - Co)k°‘^katp v - AgH ly 

+ iH- Co),ctk‘^^dp + g^P{iH - Co)V2)_^9„ + (H - Co)V2Hi, 

= -V {AgH + HiH^/2 - 2K) -b Co(2F - HCo/2)) 

+ iH- Co)ig^PH^t.dp - k^^.^dp) 

= -V {AgH + HiH^/2 - 2K) -b Co(2F - HCq/2)) . 

Here, we used the Codazzi-Mainardi equations for the second identity. 
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Appendix B. Covariant differentiation and curvature 

Let Cij be a Riemannian metric on a sufficiently smooth manifold M, and let 
denote its matrix inverse. Then the corresponding Christoffel symbols are given by 

'Tjy- = + djCu - didj), 

Furthermore, for the total covariant derivatives of a vector field Y and a (1, l)-tensor 
field T we have 

while their divergences are given by contracting i and j (or k using e). Furthermore, 
the e-symmetric part D’^Y of V^Y can be written as 

{D^Y)l = le^\eM{V^Y)\+eu{V^Y)^u)- 

Now, let e be a second metric with corresponding Christoffel symbols . Recall 
that the difference := — ^Ff^- of two connections is a tensor, and note tha10 

E = 6*^6, 

V"E = e*{V^fe+{V^ef. 

For a scalar function / we have 

(gradg fY = = (grad, ff + (?^ - 

divg Y = diY^ + = div, Y + = div, F + g * V"e * F, 

Aef = divg gradg / = A,/ + dive(gradg - grad,)/ + E)fc(gradg ff 
= Aef + {e — e) * r(e, e) * (V®)^/ + r(e, e) * V®e * V/. 
Furthermore, we have 

{D^Y)i = lg^fc(gfcKV«F)' + eA/(V^F)') + (g^ * E * F)^, 

and thus 

D®F = D®F + (g — e) * r(g, e) * V®F + r(g, e) * V®g + F 
Similarly, we compute 

divg T = div, T + (g - e) * r(g, e) * V^T + r(g, e) * V®g * T. 
Furthermore, we can prove by induction that for any tensor field T we have 

k 

(66) (V'^)'=r = {V^ fT + XI n * {{vye)p^ {V^Yt, 

Pi ,l i—1 

where pi^l G N>o and I = k. Using Riemannian normal coordinates 

with respect to e, we may assume in the following calculations that the Christoffel 


^Here and in the following, we have a slight ambiguity of notation: It is not specified wether the 
covariant tensor eij or the contravariant tensor enters the formulae. This ambiguity, however, 
is irrelevant in the analysis. 
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symbols (but not their derivatives) vanish. Hence, we compute 




52 yfc + r' + diY^ + 

-%(5,y'= + 'rLy'") 

{{V^fY)% + Y’- + d,Y^ 

+ ^rUdjY"^ + ^r^Y^} - %{diY^ + ^'rLy“) 

((v«)2y)^^. + {v^Y)%Y^ + s,^,(v'=y)( + sL((v'=y)7 + yjy^) 
- s(/(^vy)f + ELy™). 


In particular, we have 

Agy = A^y + (e - e) * (v^)^y + e * (V^E * y + E * V^^y + E * E * y) 

= Agy + (e — e) * r(e, e) * (V®)^y + r(e, e) * (V®)^e * y 
+ r(g, e) * (V®e)2 * y + r(e, e) * V®e * V'^y 

Now, let r be an orientable submanifold of M of codimension 1, and let Ve and 
r'e be equally oriented unit normal fields on T with respect to e and e, respectively. 
Employing a Gram-Schmidt orthonormalization it is not hard to see that 

i^s = Ve + {e - e) * r(e, e). 


Thus, we have 

+ E * r'g = V®r'e + ?'(e, e) * V®e. 

In view of = —eij(V® z/g)*9^, where greek and latin indices refer to coordi¬ 

nates on r and in M, respectively, we deduce 

fcg = fee + (e - e) * fee + r{e, e) * V®e, 

He = He + ie-e)* r{e, e) * + r(g, e) * V^e, 

ATg = Aet{g°‘^{ke)sp) = det {g°‘^{ke)si 3 + (e- e)* r(e, e) * ke + r{e, e) * V®e) 

= Ke + r(e, e) * ((e - e) * fc * Ve, (Ve)^). 

Finally, we consider the special case of T being a closed surface in R^, and we 
let be the graph of a height function h on T as in Subsection 13.21 Note that 

(67) (p~f^^{dAh) = \detd(ph\dA and {(ph\^)~^{dsh) = \ detd{(fhU)\ds, 

where ds and dsh denote the line elements of arbitrary curves 7 C T and fhij) C 
r?i, respectively; the determinants must be taken with respect to orthonormal bases 
in the respective tangent spaces. If we define e as in the proof of Theorem 13.131 
then we have 


( 68 ) 


\ks\e = \kh\o‘fih, He = HhO(ph, etc. 
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From (|67|) . (l68l) . and (1551) we easily obtain for p G [1, oo] 

{dAh)\ < rh\dA\ < rh\(pj^^idAh)\, 
\i‘PhU)~^{dsh)\ < rh\ds\ < rh\iph\'i)~^idsh)\, 

||fc/t||p ^ T;i(l + |/l|2,p), 

(b9) 

\\Hh\\p < rh{l + \h\2,p), 

\\^Hh\\p<rt{l + \h\3,p), 

\\Kh\\p < rh{l + \h\l2p)- 

For the fifth estimate we used the interpolation inequality dni). 


Appendix C. Function spaces 


In this section we give the definition and some properties of the function spaces 
we are using in the present paper, namely the Sobolev-Slobodetskij spaces Wp and 
the Besov spaces Bp Let e be a Riemannian metric on a smooth d-dimensional 
manifold M. For measurable tensor fields T and 1 < p < c» let 

WTVp--j \T\ldVe and ||T||oo := esssup^ |F|e, 

M 


where dVe is the volume element corresponding to e and | • |e is the norm induced 
by e on all tensor bundles. We write T G Lp(M), 1 < p < oo, if HFlIp is finite. For 
k G No and 1 < p < c» we write T G Wp{M) if (V®)*T G LP{M) for all I G No with 
I < k, and we define 


|TU,p := ||(V«)'=T||p. 

Furthermore, for s = k S with k G No, S G (0,1) and 1 < p < oo we write 
T e W^(M) if Tg W^(M) and 

KV-mx) - (V^mp)lP 


\T\ip ■■= 


M M 


de{x,y)‘i+^P 


■ dVe{x) dVeiy) < 00, 


where de denotes the distance function corresponding to e. We will use the notation 
H^{M) := IF|(M). For arbitrary real s > 0 we define 


\\T\\s,p:= E \T\i.p+\T\s,p. 

0<l<ls] 

where [s] is the largest integer smaller than s; note that this summation doesn’t 
make sense from the point of view of physical dimensions. Now, let T be a measur¬ 
able tensor field defined in M = and consider the operator ShT := T{ - -\-h) — T. 
For k gN and 1 < p, g < oo we write T G Bp^^iW’') if T G and 

\T\k,p,, ■■=\\\h\-^\\Sl^'^-^T\\p\\ 

where the Lg-norm is taken with respect to the measure dh/\h\'^. For s = k 6 
with k G No, S G (0,1) and 1 < p, g < oo we write T G R® if T G Wp(R‘^) 

and 

\T\s,p,, ■■= |||/ir'l|4V'=T||p|| ^ 


< oo. 
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For arbitrary real s > 0 we define 

Ill’ll S,P,(J \T\l^p + \T\s^p^q. 

0 </<[s] 

By the theorem of Section 2.5.12 in [45] this definition of the spaces is 

equivalent to the usual definition based on the Fourier transform and Littlewood- 
Paley theory; the latter can be extended to arbitrary s G M. For non-integer s and 
1 < p < oo we have Bpp{U.‘^) = lFp(IR'^) algebraically and topologically, since the 
semi-norms | • |s,p and | • |s,p,p coincide in this case. By Proposition 2 (iii) of Section 
2.3.2 and the theorem of Section 2.5.6 in we have even for 

integer k. Now, let M be a smooth, closed manifold with some finite atlas ipj and 
a subordinate partition of unity (pj. For real s > 0, 1 <p,q < 00 , and measurable 
tensor fields T defined on M we write T G Bp^^{M) if (pj T) o G Bp g(R‘^) for 
all j, and we define 

\\T\\s,p,q l|s,p,9, \T\s,p,q ■='^^\iPjT) O-lpj |s,p,g- 

j 3 

It is not hard to see that the definition of the spaces W^{M) by this localization 
procedure is equivalent to the one given above; hence, we have Bpp{M) = Wp{M) 
for non-integer s and 1 < p < 00 and i ?2 2 (-^) = H^{M) for integer k. Next, let M 
be a smooth bounded domain in For real s > 0, 1 < p, g < 00 , and measurable 
tensor fields T defined in M we write T G B^ q{M) if T is the restriction of a tensor 
field T G Bp^q(W^). We define 

||T|U,p,, :=inf||f|U,p,„ 

T 

where the infimum is taken over all extensions T of T. By the theorem of Section 
3.3.4 in |45j there exists a common linear and continuous extension operator E : 
Bp^g{M) —>• Bp^g{R‘^) for all 1 < p, q < 00 and uniformly bounded real s > 0. It 
follows from the theorem of Section 5.2.2 in |?5] that the definition of the spaces 
Wp{M) by this restriction approach is equivalent to the one given above; hence, 
again, we have Bpp{M) = Wp{M) for non-integer s and I < p < 00 and B 2 ^ 2 (-^) = 
H^{M) for integer k. We define 

\T\s,p,q '■= \ET\s^p^q 

for some fixed instance of E. Furthermore, by Proposition 2 of Section 2.3.3 in [5^ 
and by the construction of the spaces given above we have 

(70) Bl^^{M) ^ Bl^^{M) 

if <71 < <72 for M = and for M a smooth, closed manifold or a smooth bounded 
domain. 

Lemma C.l. Let T be a closed surface contained in a smooth bounded domain 
fl C R^. For s > 0, h G 7 > 0, and 2 < p < 00 such that a := 

{t — 2/p — l)/s G (0,1) we have 

(71) \h\t,pp<c{\h\^+2jh\\lX + \\hh,2). 

For s > 0, q G i/®(r), 0 < 7 < s, and a = tfs we have 

(72) Mq 2 <c(|q|“ 2 l|g||^-“ + lkll 2 ). 
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For s > 0, u € Fl^ifl), v = Pr[u]r G t > 0, and 2 < p < oo such that 

a := (t — 2/p + l/ 2 )/(s -|- 1 / 2 ) G ( 0 , 1 ) we have 

(73) |T|t,p ,2 < c(|i;|“+i, 2 l|Vu||^-“ + llVulb). 

For s > 0, a £ FI^{n), t > 0, and 2 < p < oo such that a := (i—3/p-|-3/2)/s £ (0,1) 
we have 


(74) |akp, 2 <c(|a|“ 2 l|a||^-“ + |H| 2 ). 

For s > 0, u G -ffg (11) n t > 0, and 2 < p < oo such that a := {t — 2/p + 

l/ 2 )/s £ ( 0 , 1 ) we have 

(75) IMrlcp .2 < c(k|“+i. 2 l|Vu|k-“ + llVulh). 


Proof. Let us begin with the proof of dzlD. From the theorem of Section 2.4.2 in 
pS] we have 

where Sq, = as -I- 2. Hence, using the embedding iL®“(R^) from the 

theorem in Section 2.7.1 in |45j with t,p as in the assertion, we obtain 

lklkp,2 < C ||/l||s^2,2lkll2,2 ■ 


In view of the definition of the Besov norms on F the same estimate is true for all 
h G iJ''+2(r); in fact, we even have CD. Indeed, if this estimate was false, there 
would exist a sequence (/i„) C i7®''‘^(r) such that 


\h 


■nls+2,2 


\\hr, 


11 —a 

l2,2 


|h„||2.2 < ^ Ikn|kp,2 < ^ lkn||“+2,2lkr, 


I 1 —Q; 

l2,2 ■ 


Dividing by ||ft ,„||2 2 “ and taking the a-th root we obtain 


|ki|s-|-2,2 + Ikn||2,2 < ) \\hn\\ s-1-2,2 ■ 

We may assume without restriction that ||/in||s-i-i ,2 = 1- Hence, we have hn —>■ 
0 strongly in i7^(r), weakly in i7®+^(r), and thus strongly in by fbe 

compact embedding i7®+^(r) ^ which follows from Remark 1 of Section 

4.3.2 in [45]; here, [s + 2] is the largest integer smaller than s + 2. Since also 
|hn|s-i- 2,2 —0 this is a contradiction to ||/in||s-i- 2,2 = 1- This proves (FfTl) . The proof 
of (|72ll proceeds essentially along the same lines. 

The proof of (I73p is very similar using 

(i7«+kK"),i7^(K'))l-a.2 = i7*“(R') S‘_ 2 (K") 

with Sa = a;(s + 1) + (1 — a)l /2 and t,p as in the assertion, the continuity of 
the trace operator u 1 —>■ [ttjr, —>■ which follows from the theorem of 

Section 3.3.3 in [45], and Poincare’s inequality. 

For the proof of ([74| we note that 

(7JkR"),T2(R"))l-a.2 = 77'*“ (K") S‘,2(R"), 


with Set = as and t,p as in the assertion, and we make use of the common extension 
operator E : i7®(f2) —>■ F[^{M?), E : 72 ( 11 ) —t 72 (R.^) introduced above. 

Finally, for the proof of (1751) we proceed very similarly by making use of 

(i7"+i(R^),i7i(R^))i-ct.2 = i7'*“(K^) .b‘;^2"(R'^). 
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where Sa = as + 1 and t,p as in the assertion, the trace operator u i—>■ [u]r, 
-^p, 2 (r), see again the theorem of Section 3.3.3 in and Poincare’s 
inequality. □ 

In the proof of Theorem 13.131 we are considering perturbations of a some fixed 
metric e. Whenever we want to emphasize the dependence of the semi-norms on 
the metric we write e\T\s,p- Now, it is not hard to see that for s > 0 and 1 < p < oo 
we have 

(76) e|T|yp < r(||e||oo, ||e-'||oo, ||e||oo, \\e-^\\oo) e\T\s.p, 

where the L°“-norms are taken with respect to a fixed background metric. 


Lemma C.2. Let M be a smoothly bounded domain or a smooth, closed manifold. 
Then, for suffieiently regular tensor fields Ti,... ,Tfc defined in M, and 6 € (0,1) 
we have 


nr. 


5,2 


i=i ‘=1 

1^3 


Tillpij \Tj\d,pjj,2 


where ptj G [ 2 ,oo] sueh that ^/Pij 


1/2 for all j. 


Proof. It is sufficient to prove the estimate in the case M = Mr. Indeed, then we 
deduce that for a closed manifold M 

k k 


nr|„<i:|(wnr)=fr‘ 


(5,2 

2=1 I 


2=1 
k k 


5,2 


^EEn WpijliPlTj) O \s,pjj,2 


I i=i 

k k 






while for a bounded domain M we have 

k k k k 


2=1 
k k 


5,2 




— E/ n 


4 = 1 *=1 

1^3 


where E is the extension operator introduced above, and the first estimate is a 
trivial consequence of the definition of the | • semi-norm. Now, let M = M'^. We 
assume for simplicity that k = 2; the general case follows by induction. We have 

|TiT 2 |y 2 = |||/ir1<5,,(rir2)||2|L an 

’NF 

< |||/l|-iTi^;,T2||2|L .. +\\\h\-^\\6hTiT2{- + h)hL an 

< \\Ti\\pJ\h\-^\\ShT2\\p,,\\^ ^ + \\T2\U\\\h\-^\\6nT,\\p,,\\^ ^ 

’FF 

= \\Tl\\q,2\T2\s,q22,2 + 11721 | 921 15 . 911.2 ■ 
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Here, we used the triangle inequality for the first estimate and Holder’s inequality 
for the second estimate. □ 

Corollary C.3. Let M be a smoothly bounded domain or a smooth, closed manifold. 

Then, for a sufficiently regular tensor field T of class {k, 1) defined in M, an analytic 

bundle homomorphism f mapping the tensor bundle of class (k, 1) to some other 

tensor bundle, and 6 € (0,1) we have 

\f{T)\5,2<f{\\T\U\TUa, 

where f is an increasing analytic function depending on f. 
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